Quantum signatures of 
Solar System dynamics 



Arkady L. Kholodenko 



Abstract Let uj{i) be period of rotation of the i-th planet around the Sun 
(or ujj{i) be period of rotation of j-th satelhte around the i-th. planet). From 
empirical observations it is known that within margins of experimental errors 
^jniCi;(j) = (or X)j = 0) fo'^ some integers rii (or rij ), different for 

different satellite systems. These conditions, known as resonance conditions, 
make uses of theories such as KAM difficult to implement. The resonances in 
Solar System are similar to those encountered in old quantum mechanics where 
applications of methods of celestial mechanics to atomic and molecular physics 
were highly successful. With such a successes, the birth of new quantum me- 
chanics is difficult to understand. In short, the rationale for its birth lies in 
simplicity with which the same type of calculations can be done using meth- 
ods of quantum mechanics capable of taking care of resonances. The solution 
of quantization puzzle was found by Heisenberg. In this paper new uses of 
Heisenberg's ideas are found. When superimposed with the equivalence princi- 
ple of general relativity, they lead to quantum mechanical treatment of observed 
resonances in the Solar System. To test correctness of theoretical predictions 
the number of allowed stable orbits for planets and for equatorial stable orbits 
of satellites of heavy planets is calculated resulting in good agreement with ob- 
servational data. In addition, the paper briefly discusses quantum mechanical 
nature of rings of heavy planets and potential usefulness of the obtained results 
for cosmology. 
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1 Introduction 



1.1 General comments 

The role of celestial mechanics in development of modern quantum mechanics 
is well described in lecture notes by Born [1] . Surprisingly, usefulness of the 
atomic mechanics to problems of celestial mechanics has been recognized only 
very recently [2,3]. Closely related to these papers is the paper by Convay at 
al [4] where methods of optimal control and genetic algorithms were used for 
mission planning problems^. In this work we extend the emerging reverse trend. 
For this purpose it was nesessary to critically reanalyzed the logical steps leading 
from the old to new quantum mechanics in the light of available astronomical 
observational data. For the sake of uninterrupted reading, some not widely 
known facts from history of quantum mechanics are presented in nontraditional 
setting involving the latest results from mathematical and atomic physics. These 
facts are very helful for formulation of the problems to be solved in this paper. 
Thus, below, we discuss some historical background information first. 

1.2 Resonances in old atomic mechanics 

In 1923-24 academic year in Gottingen Max Born replaced planned two- semester 
lecture course in celestial mechanics by the course in atomic mechanics. Con- 
trary to the standard superficial descriptions of "old" quantum mechanics which 
can be found at the beginning of any textbook on quantum mechanics, the 
achievements of " old" quantum mechanics go far beyond calculation of spectra 
of Hydrogen atom. In fact, the optical and X-ray spectra of almost all known at 
that time elements were found accounting even for the fine structure relativistic 
effects. The theory of quantum angular momenta was developed and used in the 
theory of polyatomic molecules. The effects of Zeemann and Stark were consid- 
ered as well, etc. If one would make an itemized list of problems considered in 
"old" quantum theory and would compare it with that for "new" quantum the- 
ory, surprisingly, one would not be able to find an item which was not treated 
within the "old" formalism. With such an impressive list of acomplishments 
it is hard to understand why this formalism was abruptly abandoned in favour 
of "new" quantum mechanics in 1925. To explain this, we would like to bring 
some excerpts from the paper by Pauli and Born [5] . Being thoroughly familiar 
with works by Poincare' on celestial mechanics, they were trying to apply these 
methods to multielectron atoms. For this purpose they were using methods of 
theory of perturbations to account for electron-electron interactions. By doing 
so they obtained the same types of divergencies as were known already from 
calculations of planetary dynamics. By realizing the asymptotic nature of the 
obtained expressions, they decided that to keep just few terms in these expan- 
sions is the best way to proceed. By doing so a reasonably good agreement with 
experimentally known location of spectral lines was expected to be obtained. 
Such a state of affairs had caused frustration for Bohr who conceded that only 

^That is problems involving optimal interplanetary travel in Solar System. 
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those dynamical systems which admit a complete separation of variables are 
quantizable^ . If such a separation is absent, according to Bohr's current opin- 
ion, the system should not possess a discrete spectrum so that visible lines in 
spectra of elements other than Hydrogen should/must be much wider. On the 
theoretical side such an assumption calls for development of methods enabling 
to determine the widths of spectral lines and of distribution of the intensity 
within these widths. Such an intensity is expected to be connected with the 
underlying mechanical motion inside the atomic system. 

Spectroscopical data for almost entire periodic system were readily available 
at the turn of the 20th century [6] . Bohr was well aware of these data and used 
them for his search for correct atomic model (along with Rutherford's results 
of 1912 on scattering from the Hydrogen atom). In particular, he looked at the 
data for Helium in 1913 and published his findings in Nature [7]. For the sake 
of arguments which will follow, we found it helpful to reproduce some of the 
data from his Table 1 below. 



Table 1 
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These data were compared with those for the Hydrogen for which he used the 
analogous table (Table 2)"'. For some reason, the data in his Table 2 did not 
contain the error column. Since the wavelength A in both cases was measurable, 
it was possible to evaluate the ratio K/j/K/fe, where K= A(^ — ^j-) • 10^°, which 
was found to be 4.0016. At the same time, Bohr's own calculations gave for K 
the following value: K= 2^2^t^]^^ , with h being the Planck's constant, Z and 
M being the charge and the mass of the nucleus, c being the spc;ed of light 
and e and m are being the charge and the mass of the electron. By assuming 
Mffe = 4Mjy and Zhb = 2Zjy , one readily obtains for Kjy/Kjje the result: 
4.00163. It is in good agreement with that obtained experimentally. In doing 
such calculations Bohr assumed that each electron in Helium can be treated 
as if it is a Hydrogen-like. This surely implies that the width of spectral lines 
for Helium should be practically the same as those for the Hydrogen atom. 
Nowadays we know [8] that all atomic spectra have some finite linewidth. This 
linewidth is determined by factors such as: a) the coUisional broadening, b) the 
Doppler broadening and c) the natural broadening. Each of these is having some 

^E.g. read [5]. 
Which we do not reproduce. 
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further ramifications. Hence, from the standpoint of modern knowledge one can 
interpret Bohr's conclusions made in 1922 as acknowledgement of the fact that 
spectra of elements other than Hydrogen are broader because of natural reasons 
(so that one should take into consideration the data from the error column in 
Table 1) without invalidation of Bohr's major quantization assumptions. It 
should be noted though that such a conclusion leads to the question: why the 
very same factors are affecting the Hydrogen atom much less? 

Unhappy with his conclusions, Bohr asked Born and Heisenberg to make 
more rigorous calculations for Helium using perturbational methods analogous 
to those developed in the paper by Pauli and Born in 1922. Their findings were 
published in 1923 and resulted in practically total failure in accurate determi- 
nation of energies of the ground and excited states for Helium atom. This fact 
is documented in Bern's lecture notes [1]. 

Since the Helium atom calculations were made by Heisenberg (under Horn's 
supervision) it might be not too surprising that, after all, it was Heisenberg 
who found the way out of the existing difficulties. The logic of his reasonings 
is discussed from the modern mathematical point of view in Section 2. In 
this introductory section we would like only to put his work in some historical 
perspective. For this, we need to make few comments on his joint work with 
Born. For the sake of space, we refer our readers to the cited literature for 
details. 

1.2.1 3-body problem and the He atom (old quantum mechanics 
results) 

The unperturbed Hamiltonian for He was chosen as H =-A(J]^^+J2^^) with 
the constant A — 2n'^e^m,Z'^ while the perturbation was chosen as Hi = e^/R 
with Ji and J2 being the Bohr-Sommerfeld (B-S) adiabatic action integrals for 
electron 1 and 2 and R being the Euclidean distance between them. After H, 
He is the first nontrivial 3-body mechanical system whose behavior is amended 
in accordance with the rules of old quantum mechanics by requiring both of 
these integrals (i.e. Ji and J2 ) to have their lowest value, i.e. h, so that the 
unperturbed energy for He is twice that for H. Since the energy W for H is known 
to be = — the frequency w of rotation of the electron at its stationary 
orbit in the action-angle formulation of classical mechanics is obtained as w = 
^ In obtaining this result it was assumed that J is continuous 

variable and, only after the differentiation is done, J is assumed to be discrete: 
J= nh. It is important to realize at this point that exactly the same logic was 
used in Heisenberg's paper on quantum mechanics to be discussed in Section 
2^. Hence, for Helium atom within the approximations made the rotation 
frequencies of both electrons are the same. This fact is known in mechanics 
literature as accidental degeneration. In view of its crucial importance for this 
paper, we would like to pause now in order to provide more accurate definitions. 

"'E.g. see page 140 of Ref.[l]. 
^E.g.see equation.(19) below. 
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In terms of the action -angle (!,</?) variables the Hamilton's equations of 
motion for a completely integrable system can be written as: 




(1) 



where the boldface indicates that the dynamical system with Hamiltonian H is 
multicomponent (in general case). Solutions of the system (1) are: I, =Ci,Lpi = 
LOi{l)t-\-Ci , i = 1~ N. It is assumed that Cj and Cj are some known constants. 
In view of this result, any mechanical observable ^(p, q) made of generalized 
momenta p and generalized coordinate q can be Fourier decomposed as 



where n = {m, ...,nAr}. Accordingly, n • ip =X]i=i'^j'/'i- Such a Fourier decom- 
position is expected to exist even for those perturbed systems for which the 
empirically observed orbits are closed (as seen in the case of planets in our So- 
lar System). Dynamical system is considered to be accidentally degenerate if 
the relation 



holds for some fixed set of integers n and is degenerate if (3) holds for any set 
of n's. In this work this condition will be alternatively called as resonance 
condition in accord with modern terminology. 

1.2.2 3-body problem and the He atom (modern perspective) 

Failure of methods of celestial mechanics (modified by the B-S quantization rule) 
to accurately compute the ground and excited states for He has led Heisenberg 
to discovery of his matrix quantum mechanics in 1925. The problem appeared 
to be completely solved until 1960 when some difficulties emerged when the 
standard Hartree-Fock type variational calculations become inadequate for de- 
scription of doubly excited electron states [9]. The same authors notice that 
by 1990 the improvements which were made in 60ies failed again, especially for 
the extreme excitation regime which cannot be described using single electron 
quantum numbers. The way out of existing difficulties was associated with ac- 
curately designed semiclassical methods. The backbone of these semiclassical 
descriptions "are the periodic orbits of the full classical two-electron system 
without any approximations" . Tanner at all [9] notice (e.g. read page 523) that 
"The classical two-electron atom is neither integrable nor fully chaotic. The 
apparently regular spectrum as well as the breakdown of approximate quantum 
numbers for highly doubly excited states and the enormous variation in the de- 
cay widths for resonances can be understood by studyng classical mechanics in 
detail. Qualitative results can be obtained by exploiting semiclassical periodic 
orbit theory." In other words, use of classical mechanics is quite sufficient for 
determination of both the ground and excited states of He and only for extreme 



oo 



^ = J2 ^nexp(m-^), 



(2) 



(3) 
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case of highly doubly excited He the description becomes qualitative. Classi- 
cal dynamics of He used for calculations of spectra is essentially the dynamics 
of the restricted 3-body problem^ and, as such, exhibits chaotic and regular 
regimes. Since the experimental information which can be deduced from the 
low lying spectral excitations of He does not allow to disentangle regular and 
chaotic parts of dynamics, it has become possible to simplify things further by 
recalculating spectra of He, Li, Be. and diatomic molecules made out of these 
and other atoms, and also of H2, by revisiting Bohr's 1913 calculations [10-12]. 
These calculations involve a simple minded minimization of classical functionals 
of the type considered by Bohr in 1913. Mathematical justification of such a 
procedure was found by Chen et al [13] The accuracy of results obtained with 
help of such classical calculations (employing however the B-S quantization rule 
!) compares well with incomparably more elaborate traditional quantum me- 
chanical calculations. It should be noted though that in his Nobel prize winning 
address Bohr (1923) was talking about his great success in calculating spectra 
of almost all elements of the periodic table using the B-S quantization rule and 
simple minimization procedure. It took another 80 years or so to bring these 
calculations to the level comparable with the best known quantum mechanical 
calculations! 

1.3 Resonances in celestial mechanics 

On page 265 of his lecture notes [1] Born writes: "Accidental degeneration 
is a rare and remarkable exception in astronomy; the odds against (Eq.(3)) 
being exactly fulfilled are infinite. A close approach to it is found in the case 
of perturbations of some minor planets (Achilles, Patroclus, Hector, Nestor) 
which have very nearly the same period of revolution as Jupiter. In atomic 
theory, on the other hand, where J^'s can have only discrete values, accidental 
degeneration is very common.^^ As result of such an accidental degeneracy 
Heisenberg's attempt at perturbative calculations for He failed miserably. Such 
a failure caused him to reconsider the whole computational scheme resulting in 
an ultimate breakthrough in 1925 leading to new quantum mechanics. 

Before discussing his contributions from the modern perspective, we would 
like to make few remarks regarding the accuracy of astronomical data in Born's 
lectures. In 1968 Molchanov [14], while analyzing the astronomical data, came 
to the conclusion that the accidental degeneracy for Solar (and, very recently. 
Solar-like [15,16] system(s) is as common as in atomic systems. In Table 2 
(below) taken from his work (Molchanov's Table 1) we reproduce some data 
taken from this rolbroiico. 

^Modified by the fact that in the atomic case electrons repel each other and have the same 
mass. 
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Table 2 





Planet 








ni 


n2 


na 


n4 


n5 


ne 


ny 


ns 


ng 


1 


Mercury 


49.22 


49.20 


0.0004 


1 


-1 


-2 


-1 

















2 


Venus 


19.29 


19.26 


0.0015 





1 





-3 





-1 











3 


Earth 


11.862 


11.828 


0.0031 








1 


-2 


1 


-1 


1 








4 


Mars 


6.306 


6.287 


0.0031 











1 


-6 





-2 








5 


Jupiter 


1.000 


1.000 


0.000 

















2 


-5 








6 


Saturn 


0.4027 


0.400 


0.0068 














1 





-7 








7 


Uranus 


0.14119 


0.14286 


-0.0118 




















1 


-2 





8 


Neptune 


0.07197 


0.07143 


0.0075 




















1 





-3 


9 


Pluto 


0.04750 


0.04762 


-0.0025 

















1 





-5 


1 



For satellite systems of Jupiter, Saturn and Uranus Molchanov's paper also 
contains tables similar to our Tabic 2. According to the book by Bclctsky [17], 
in view of the resonance nature of our Solar System, uses of KAM theory [18] 
for explanation of planetary stability typically fail. 

To understand the data in Table 2 several comments arc in order. First, the 
displayed frequencies are measured in the system of units in which the Jupiter's 
frequency was chosen as the unit of measurement. Second, in view of Eq.(3), the 
first row of data from Table 2 should be actually read as wi — a;2 — 'i.uJ?, — (x'4 = 0. 
All other rows should be treated accordingly. The theoretical frequencies ujj 
are those which satisfy the resonance conditions exactly while wf denote the 
observed frequencies. The data for Pluto should not to be considered in terms 
of resonances for the following reason. 

Consider a scalar product n • = A in Eq.(2). This can be looked upon as 
representation of the vector A in the coordinate basis {95}. The coordinate basis 
can be changed with help of some matrix A so that A = n-A-<^. It can be argued 
[1],[14] that detA = 1, so that the matrix A must be a unimodular square 
matrix. Only for the sake of this requirement the data for Pluto in Table 2 were 
assigned in a way given in the Table 2. Next, Aco /lu should be understood as 
{u)^ — uf)/u^. After this, the obtained error margins can be compared against 
those for He in Table 1. Such a comparison indicates that the accuracy in both 
cases is essentially the same. It is such that Bohr was able to obtain using his old 
quantum mechanical theory a reasonably accurate ratio /Kh^ in agreement 
with experiment. It makes physical sense to blame the intrinsic inaccuracy of 
the collected data (e.g. that in Table 1) for observed frequency discrepancies. 
Hence, along with Bohr, it is reasonable to claim that, with exception of H, 
other atomic systems are not quantizablc bccaiisc of these discrepancies . The 
same reasoning should then be applied to the planetary systems, especially in 
view of critique of Molchanov's work by Henon [19] and Backus [20]. These 

^Subsequent developments of quantum mechanics demonstrated that Bohr was apparently 
wrong. We say "apparently" in view of the results of Section 2 where correct quantization 
prescription is discussed based on improvement of Heisenberg's ideas. In view of results of 
Section 2, it is reasonable to say that Bohr's intuition was nevertheless correct but the 
situation can be improved rigorously using mathematical methods which were not available 
in Bohr's time. 
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authors argued that the error margins in Molchanov's tables arc too large for 
the resonances to be considered seriously. The comparison between the Tables 
1 and 2 indicates that even though the arguments by Henon and Backus may 
be mathematically correct, they do not have sound physical support due to 
intrinsic inaccuracies in measurements which cannot be substantially improved*. 
Thus, Molchanov's data and their interpretation remain correct at the physical 
level of rigor even without additional explanations made by Molchanov [21-22] 
in defence of his results. Furthermore, subsequently obtained results by Brin 
[23] and Patterson [24] strongly suggest that, at least pairwise, planets are in 
resonance with each other. This is true in particular for the heavy planets in 
our Solar System as demonstrated by Ferraz-Mello et al [16]. Evidently, the 
linear combination of such pairwise resonances leads back to the Molchanov- 
type results. 

It should be noted though that quantization prescriptions discovered by 
Heisenberg remain correct even in the case when there are no resonances: They 
rely only on the existence of the closed stable orbits. Existence of resonances, in 
fact, simpliGes matters considerably since it makes the task of establishing the 
quantum-classical correspondence much easier as explained in the rest of this 
paper. 

Prior to Molchanov's 1968 work an effort to explain the ubiquity of reso- 
nances in Solar System using methods of classical mechanics was made by Gol- 
dreich [25] (1965) who demonstrated that "special cases of commensurate mean 
motions are not disrupted by tidal forces" . Moreover, he proposed that it is the 
tidal forces which drive otherwise incommensurate system to commensurability. 

Thus, the problem of stability of our Solar System is very much the same as 
that for the multielectron atoms. In both cases the accidental degeneracies (res- 
onances) preclude systematic use of standard perturbational methods. Unlike 
more traditional classical mechanics treatments [26] , we apply Heisenberg-style 
arguments ultimately aimed at explanation of Solar System stability. For the 
sake of space, we are not discussing in this work the spin-orbit- type resonances 
also ubiquitous in the Solar System [26]. 

Finally, we would like to mention that development of our quantum me- 
chanical formalism proceeds in historical accord with that for atomic systems 
for which the static (spectral) problems were considered first. The dynami- 
cal problems of atoms/molecules formation as well as their stability towards 
disintegration were considered only afterwards. Hence, only the spectral-type 
problems will be discussed in this work. 

1.4 Organization of the rest of this paper 

Existence of stable closed orbits, of resonances, as well as the lack of dissipa- 
tion (in spite of presence of tidal effects) in Solar (and Solar-like) System(s) are 
indicative of quantum nature of the orbital motions in the Solar System. Nev- 
ertheless, the formalism of quantum mechanics in its traditional form present 

^E.g.read Sections 3 and 4 where these measurements are discussed for objects such as 
Solar System, etc. 
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in textbooks for students cannot be used. To apply methods of quantum me- 
chanics to celestial mechanics is possible with use of Heisenberg's original ideas 
updated with help of the latest mathematical results. Sections 2 and 3 as well 
as Appendices A and B provide a self-contained overview of quantum mechanics 
based on Heisenberg's ideas. They provide needed background for the actual 
quantum calculations in celestial mechanics which are performed in Section 4 
(supplemented with Appendix C). The main results of this section (and the 
whole paper) are summarized in Table 4. In this table the number of stable 
orbits for planets of Solar System as well as the number of stable orbits for 
satellites of heavy planets (Jupiter, Saturn, Uranus and Neptune) is calculated 
and compared against the observed numbers. Unusually good agreement be- 
tween the calculated and observed numbers for Solar System, and the satellite 
systems of Jiipitcr and Saturn is obtained resulting in further siiggestions for 
observational astronomy of the Solar System. In Table 3 a comparative sum- 
mary of main theoretical assumptions of both quantum atomic and quantum 
celestial mechanics is given. Using it, the motion of rings aroimd heavy planets 
is studied briefly resulting in the same conclusions about the quantum nature 
of such type of motion. Some auxiliary mathematical results needed for these 
calculations arc presented in Appendix D. In study of all cases, including dy- 
namics of planetary rings, it follows that the equivalence principle of general 
relativity plays the decisive role in development of quantum celestial mechanics. 
This fact caused us to write Section 5 in which the effects of general relativity 
on quantum mechanics of Solar System are studied. In it we discuss how the 
obtained results (for which the importance of the Lorentz group S0(2,l) is em- 
phasized) should be amended if we are interested in knowing to what extent 
the (larger scale) symmetries of space-time typically considered in cosmological 
models of general relativity may affect the quantum dynamics of Solar System. 
Such an information can be used in reverse for probing symmetries of space-time 
at scales comparable or larger than that for our Solar System. Finally, Section 
6 is devoted to some concluding remarks. 



2 Heisenberg's honeycombs and resonances 
2.1 General comments 

In this section we discuss Heisenberg's ground breaking paper [27] on quantum 
mechanics from perspective of modern mathematics. We begin with observa- 
tion that the Schrodinger equation cannot be reduced to something else which 
is related to our macroscopic experience. It has to be postulated.^ On the 
contrary, Heisenberg's basic equation from which all quantum mechanics can be 
recovered is directly connected with experimental data and looks almost trivial. 

^Usually used appeal to the DeBroigle wave-particle duality is of no help since the wave 
function in the Schrodinger's equation plays an auxiliary role while the De Broigle waves are 
assumed to exist in real space-time. 
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Indeed, following Bohr, Heisenberg looked at the famous equations for energy 
levels difference ^ 

Lo{n,n- a) = -{E{n) ~ E{n- a)), (4) 

where both n and n — a are some integers. He noticed that this definition leads 
to the following fundamental composition law: 

u){n — /3,n — a — (3) + u}{n, n — (3) = w{n, n — a — (3). (5a) 

Since by design u){k,n) = —Lj{n,k), the above equation can be rewritten in a 
symmetric form as 

u}{n,m) + u}{m,k) + co{k,n) = 0. (5b) 

In such a form it is known as the honeycomb equation (condition) in current 
mathematics literature [28] where it was rediscovered totally independently of 
Heisenberg's key quantum mechanical paper and, apparently, with different 
purposes in mind. Connections between mathematical results of Knutson and 
Tao and those of Heisenberg were discovered in the recent paper by Kholo- 
denko[29]. In this work some results of this paper will be used. 

In particular, we begin by noticing that Eq.(5b) due to its purely combi- 
natorial origin does not contain the Plank's constant h. Such fact is of major 
importance for this work. In particular, the simplest resonance condition en- 
countered in celestial mechanics 

niUJi + n2UJ2 + n^LO^ — (6) 

can be equivalently rewritten in the form of (5b), where u>{n,m) = uj„ — Um- 
It would be quite unnatural to think of the Planck's constant for this case. 
Even though, the resonance condition is equivalent to Heisenberg's quantization 
condition, Eq.(5b), the reverse may not be true since frequencies in Eq.(5b) 
may be irrational. It should be noted though that such irrationality would be 
very difficult to detect experimentally in view of natural causes leading to the 
line broadening mentioned in the Introduction. Thus, from the experimental 
standpoint equations (5b) and (6) are equivalents^. Furthermore, by assuming 
irrationality wc would run into difRculty with obtaining the scmiclassical limit in 
which (it is believed) the old quantum mechanics based on the Bohr-Sommerfeld 
method of quantization should be applicable. These arguments imply that, at 
least scmiclassically, dynamics of all quantum mechanical systems is resonant. 

Equation (5b) is the basic building block of the honeycomb structure encod- 
ing all information about the spectra of quantum system. Details leading to 
construction of this combinatorial structure are summarized in the paper by 
Knutson and Tao [28] . They were used in Kholodenko's paper [29] in which some 
physical applications absent in Knutson's-Tao paper are discussed in detail. 

^"Although to maJie the frequencies independent (and, hence, irrational) is easy mathemat- 
ically, it is unrealistic to detect such fact experimentally. In other words, even though the 
critique of Molchanov's paper [14] by Henon [19] and Backus [20] could be mathematically 
justified in spite of Molchanov's counter arguments [21,22], it is only of academic value. 
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To describe such honeycomb structure in a nutsheU, let us choose the basic 
Y- shaped tripod whose edges are labeled by frequencies uj{n,m)'s is such a 
way that the total sum of these labels is equal to zero, as in Eq.(5b). The 
honeycomb is made of collection of such tripods placed on a 2- dimensional 
plane and joined with each other in such a way that the frequencies at the 
edges match. Several additional rules were set up by Knutson and Tao and 
are given in their original papers [30, 31]. Our readers encouraged at this 
point to consult the interactive web site designed by Tao [32] in order to get 
a feeling of honeycombs as combinatorial objects. For physical applications, 
other than those discussed in this paper, our readers are referred to the paper 
by Kholodenko [29]. Provided references allow us to squeeze to the absolute 
minimum the amount of mathematical information in this paper. 

With account of these remarks, we proceed with development of Heisenberg's 
arguments. In his paper of October 7th of 1925, Dirac[33], being already aware 
of Heisenberg's key paper^^, streamlined Heisenberg's results and introduced 
notations which are in use up to this day. He noticed that the combinatorial 
law, Eq.(5a), for frequencies, when used in the Fourier expansions for compo- 
sition of observables, leads to the multiplication rule a{nm)b{mk) = ab{nk) for 
the Fourier amplitudes for these observables. In general, in accord with Hcisen- 
berg, one expects that ab{nk) ^ ba{nk). Such multiplication rule is typical for 
matrices. In the modern quantum mechanical language such matrix elements 
are written as < n | O | m > exp(ia;(n, m)t) so that Eq.(5b) is equivalent to the 
matrix statement 

<n\Oi \m >< m\ k> exp{iLj{n,m)t) exp{iLo{m, k)t) 

=< n I 6162 I k > exp(iw(n, k)t) (7) 

for some operator (observables) Oi and O2 evolving according to the rule: 
6k{t) = u6kU-^,k = 1,2, provided that U^'^ = exp(-z-|t). From here it 
follows that I m >= exp(— -^f) | to > if one identifies H with the Hamil- 
tonian operator. Clearly, upon such an identification the Schrodinger equation 
can be obtained at once as is well known [34] . We shall avoid such a pathway (at 
least at this stage), however. Moreover, we also shall avoid use of Heisenberg's 
equations of motion 

ih^d = [d,H]. (8) 

Our readers may ask at this point: why it is necessary to do so? And, if 
this is the case, what else is left from the traditional formulations of quantum 
mechanics which still can be used? The answers can be found in [29,35]. For 
the sake of uninterrupted reading they are summarized also below. 

Following Heisenberg's philosophy, we shall assume that there is a set of 
classical observables {Oi{t)} which is assumed to be complete in the sense that 
the composition of any two of these observables is given by the classical fusion 

^^This paper was sent to Dirac by Heisenberg himself. 
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rule: 



(9) 



where C^^ are some known constants and {,} represents the Poisson brackets 
of classical mechanics. Accordingly, quantum mechanically, instead of Eq.(8), 
we need to consider the decomposition 

[0,H]=Y.kC^A (10) 

valid for any t ! Under such circumstances (quantum) dynamics formally disap- 
pears! This is, of course, an exaggeration since not all systems possess needed 
symmetry so that the fusion rule, Eq.(9), may not exist^^. When it docs exist, 
such an observation can be strengthened due to the following chain of arguments. 
In mathematics, expressions like Oi{t) = UOiU~^ = Ad^Oi define an orhit for 
the operator Oi in the Lie algebra (made of operators {O^}) so that the motion 
is caused by the action of elements IJ from the Lie group associated with such 
an algebra. Following the existing rules and notations in mathematics of Lie 
groups and Lie algebras [36], wc write ad^O for [O, H]. This requires us to use 
the r.h.s. of Eq.(lO) instead of the formal symbol iti-^O used in Heisenberg's 
mechanics. Evidently, we can obtain the same (or even greater) information by 
working with Ad operators instead of ad. In particular, it is useful to consider 
the trace, i.e. tr{AdjjOi} = x(Oi), which is just the character of Oj. It is time- 
independent by design. If there is no time evolution then, superficially, nothing 
happens. This is not true, however, as was recognized long time ago by Dirac 
[34]. In Chapter 9 of his book he writes : " The Hamiltonian is symmetrical 
function of the dynamic variables and thus commutes with every permutation. 
It follows that each permutation is a constant of motion. This happens even 
if the Hamiltonian is not constant^''." At this point it is important to recall 
the famous theorem by Caley [37] which states that "every finite group is iso- 
morphic to some permutation group" . It should be noted that in mathematics 
literature the "permutation" group has the same meaning as "symmetric" group 
Sn^'^- In physics and, especially, in quantum mechanics, the symmetric group 
can be infinite dimensional. The theory of such groups was unknown to Dirac 
since it was developed only quite recently [38]. This fact explains why it have 
not been in use in the traditionally written textbooks on quantum mechanics. 
Fortunately, for the purposes of this work, it is suHicient to use only a tiny 
fraction from the theory of symmetric groups. 

^^Sincc all obscrvablcs arc made of p and q variables, such a rule does exist if we decom- 
pose these obscrvablcs into power series in p and q. In those cases when such a series is 
infinite, normally, one should expect loss of integrability and, hence, loss of quantization. For 
one dimensional many body systems the situation might be repairable for suitably chosen 
interaction potentials. We shall elaborate on this remark further below, in Section 2.5. 

^^I.e. time-dependent. 

i*Here n denotes the number of elements in the group. 
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2.2 Some useful facts about S. 



In view of the fact that the character x{Oi) = tr{Ad^Oi} is manifestly time- 
independent, the orbit Ad^yOi is caused by permutations^^. These can be an- 
alyzed using methods of algebraic geometry [39] and theory of linear algebraic 
groups [40]. A brief and self-contained introduction to these topics can be found 
in Kholodenko [35] . The key concept in this field is the notion of the torus action 
T. It is directly connected with the notion of the Weyl-Coxeter reflection group 
W = N/T in which the numerator A'^ refers to some permutation group and the 
denominator T refers to those group elements (fixed points) which remain un- 
affected by permutations. Representations of Lie algebras (including the afFine 
Lie algebras) associated with these Weyl-Coxeter reflection groups produce all 
Lie algebras known in quantum mechanics and in conformal field theories [41]. 
To simplify matters, we choose another pathway in this work to arrive at the 
same results. It is better adapted for connecting the experimental data with 
theoretical constructions. 

We begin with observation that the representation theory for Sn can be 
built using representation theory for general linear group GL(N,C) acting in 
the complex space made of n copies of C^, i.e. • • -^C^ . This fact 

is known as the Schur-Weyl duality [42]. The Schur functions (to be defined 
below) are characters of GL(N,C). They play the key role in developing the 
representation theory of Sn in which both N and n can become infinite. 

Next, we recall that a partition A (finite or infinite) is a sequence 

A = {Ai,A2,A3,...} (11) 

made of integers^^ such that Ai > A2 > A3 > • • • > 0. The weight of A is denoted 
by |A| = X^^Aj. If X's are integers, and if |A| = n we say that A is a partition of n. 
Let A be a partition. It is useful to associate with it a monomial = x^^Xj^ • •• . 
Next, we introduce a symmetric function m\ as a sum of all distinct monomials 
that can be obtained from x"** by permuting of all arguments. Using these 
results it is possible to prove [43] that the Schur function sx can be represented 
with help of m\ as follows 

s\=mx+ E Kxfj,mi^. (12) 

To explain the meaning of the Kostka number Kxfi in (12) we should mention 
the one-to -one correspondence between the partitions and the Young tableaux 
[39]. In terms of such a correspondence the Kostka number Kxf^ is just the 
number of semistandard tableaux with shape A and weight /x. Hence, for not 
too large tableaux such a number can be straightforwardly computed. The 

Since, according to Dirac, the permutation operator commutes with the Hamiltonian. 
Since these numbers normaUy are identified with the eigenvalues of some matrix (finite 
or infinite) [28,29], one can relax the condition that X'^s are integers and make them rational 
or even irrational numbers but the nonnegativity and the ordering are essential. 



13 



Schur functions possess a remarkable orthogonality property. For partitions A 
and /U and properly defined scalar product <, > one can write 



(13) 



in accord with general theory of characters and, in particular, of characters of 
Sn [43]. With such defined orthogonality property of s'^s one can proceed with 
the composition (fusion) law for Schur functions. It is given by 



where |A| + |/i| = and, in view of (13), the Littlewood- Richardson (L-R) 
coefficient C^^ can be formally defined as C^^ = < sx ■ s^,Si, > .These coeffi- 
cients play an important role in representation theory of Sn analogous to the 
role the Clebsch-Gordan coefficients play in the representation theory for spin 
and angular momenta. The L-R coefficients can be obtained very easily with 
help of the honeycomb construction as discussed by Knutson and Tao [28] and 
by Kholodenko [29]. For completeness, we provide a brief sketch of how this 
can be done. 

We begin with the 1-honeycomb. It is just the Y-shaped tripod as discussed 
already. When constructing the 2-honeycomb in the plane we shall follow the 
rule that the labels for the edges of this new honeycomb should be geometrically 
and combinatorially arranged in the same way as those for the 1-honeycomb. 
This requires us to use yet another two tripods which can be joined together 
and with the third tripod only in one way in view of the imposed rules^^. Thus, 
instead of just one boundary label, e.g. Ai, in the North- West direction, now we 
shall have two, say, Ai and A2. The same applies for the South and the North- 
East directions. Thus, all larger honeycombs will have only the boundaries 
in the directions just mentioned which are labeled by the partitions A, u and 
fi. Unlike the 2- honeycomb for which the boundary labels determine such a 
honeycomb uniquely, for larger honeycombs this is no longer true. For the fixed 
set of boundary labels, normally, there will be more than one honeycomb with 
these labels. On page 1053 of Knutson and Tao paper [30] the following theorem 
is proven: Let A, and f be three pre assigned (boundary) partitions for the 
fc-honeycomb. Then the number of different honeycombs with such pre assigned 
boundary conditions is given by the L-R coefficient C^^. 

Summarizing, we have defined a set (finite or not) of mutually orthogonal 
Schur polynomials which by design forms the Hilbert space. The partitions 
and the energy levels can be put into one-to-one correspondence by using the 
honeycomb condition, Eq.(5). Such Hilbert space is designed using experimental 
data. We can look at different portions (segments) of the spectra and study their 
overlaps thanks to the composition rule, Eq.(14)^^. Unlike more traditional 
formulations of quantum mechanics requiring objects of classical mechanics as 

^^E.g. see Fig.2 in Kholodenko's paper. [29]. 

^®Very much like it is done in the case of determination of the entire DNA structure from 
its fragments. 
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an input, no reference to the objects of classical mechanics was made thus 
far. In the next (sub)sections we shall discuss the extent to which such a 
way of developing quantum mechanics is advantageous as compared with more 
traditional formulations. 

2.3 Prom combinatorics to physics 

In this subsection we follow logic of Heisenberg's paper once again. In this paper 
Heisenberg was also concerned with proper interrelation between the objects of 
classical and quantum mechanics. Naturally, he focused his attention at the 
Bohr-Sommerfeld (B-S) quantization rule 

§pdq = nh, n = 0,1,2,..., (15) 

since this rule was the only one available link between the new and old mechan- 
ics. He argued that such a rule is not exact! It is determined with accuracy up 
to a constant (unknown at the time of writing of his paper). If such a constant 
would be known, the B-S rule would become exact, that is valid for any n's. 

Prom the point of view of our present understanding of quantum mechanics 
Heisenberg's intuition was correct: the old fashioned B-S rule is valid rigorously 
only in the limit of large n's while the calculation of the constant can be done, 
for instance, with help of either the WKB or considerably more sophisticated 
theory of Maslov indices [44]. As much as these arguments are plausible, they 
are nevertheless superficial as can be found from reading the page 246 of the 
book by Arnol'd [45]. Using this reference it follows that already at the classi- 
cal level the adiabatic invariant §pdq is determined only up to some constant. 
This observation makes Heisenberg's arguments less convincing. Nevertheless, 
following Heisenberg, we assume that if the B-S quantization rule is corrected, 
it would make sense fully quantum mechanically. Presumably, under such cir- 
cumstances one can get an additional information out of it. For this purpose, 
Heisenberg introduces the Fourier decomposition of the generalized coordinate 
q as 

oo 

q{n,t) = ^ aa{n)eTCi){iijo{n,a)t), (16) 

a=— oo 

where, in anticipation of its quantum mechanical use, it is written with respect 
to some pre assigned energy level n. Using Eq.(16) the velocity can be readily 
obtained as follows 

oo 

q{n,t) = ^ iaa{n)uj(n,a)exp{iijo{n,a)t). (l^a) 

a= — oo 

The calculation of the velocity square over the total period is given therefore by 

oo 

§dt[q(n,t)f = 2-K J2 \aa{n)\ (j{n,af. (17b) 
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In view of this result, the B-S adiabatic invariant can be rewritten as 
§pdq = §mqdq = §mq^dt = 2'nm XI |aa(^)| ijjin, of = nh + const. (18) 

a= — oo 

Next, Heisenberg proceeds as follows. Since the const is unknown, it is of interest 
to obtain results which are constant-independent. At the same time, since the 
result, Eq.(18), is assumed to be exact, we have to use instead of scalars |aQ,(n)|^ 
the matrices in accord with Eq.(7). This causes us to use matrices of the type 
|a(n, n + a)|^ and |a(n, n — a)\^ depending on the actual sign of a. In addition, 
he had silently assumed that the n -dependence for amplitudes is much weaker 
than that for the frequencies Lu{n,n + a) and uj{n,n — a) so that it can be 
neglected completely. Under such conditions he treats n as continuous variable 
and differentiates both sides of Eq.(18) with respect to n thus obtaining the 
following result: 



/i = 47rTO {|a(n, n + a)| u!{n,n + a) — \a{n,n— a)\ u!{n,n — a)}. (19) 

a=0 

Obtained result takes into accoimt that Lu{mn) = —uj{nm). The validity of this 
result depends upon additional assumption about the ground state energy. If 
no represents such a state, then one must require that a(no, no — a) = for all 
a > 0. When (19) is used in combination with the results from Appendix A, 
the famous commutation rule 

[x,p]=ih (20) 

is obtained. Prom the above derivation and results of Appendix A several con- 
clusions can be drawn. 

First, the number of x-p commutators by construction is in one-to one corre- 
spondence with the number of the B-S adiabatic invariants. This means that the 
system which is completely integrable classically can be completely quantized. 
However, if classically system is nonintegrable, one cannot write the classical 
Hamiltonian and to replace x's and p's in it by the corresponding operators 
obeying commutation relations, Eq.(20), for each generalized degree of freedom. 
Although this prescription is used routinely in the existing textbooks on quan- 
tum mechanics, rigorously speaking, one cannot write down the Schrodinger's 
equation in such a case so that formally Bohr's intuition was correct. 

Second, Eq.(19) assumes that the underlying mechanical system, when it is 
written in terms of the action-angle variables, is essentially the set of indepen- 
dent harmonic oscillators. Heisenberg's derivation explicitly assumes that quan- 
tum mechanically there is a ground state-typical for the harmonic oscillator- but 
otherwise the spectrum is boundless. If the system is nonintegrable, again, the 
commutation rule, Eq.(20), is not justified. Hence, once again, one cannot write 
the Schrodinger's equation. To by pass this difficulty Heisenberg developed per- 
turbation theory (for one dimensional case only!) which uses the classical per- 
turbation theory as an input modified by the imposed (quantum) requirements 
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on amplitudes and frequencies. Although he did not discuss the resonances, 
he did checked that the obtained perturbative expressions for energies are in 
agreement with the basic equations (4) and (5) assuring correct quantization. 
No attempts to study the multidimensional case was made. 

Third, as results of Appendix A demonstrate, the experimental justification 
of the commutator rule, Eq.(20), is based on the validity of results of the first 
order perturbational calculations. Mathematically, such a procedure is ques- 
tionable or, better, might be totally unacceptable. 

Furthermore, the B-S quantization cannot be used for spin quantization 
(since formally there is no classical analog of spin, i.e. the B-S rule does not 
account for half integers). The spin has no place in the Schrodinger formalism, 
and, apparently, there is no room for spin in Heisenberg's formalism as well. 
Fortunately, this happens only apparently as we would like to discuss now. This 
is possible only because the facts just listed do not affect the main Heisenberg's 
quantization postulate, Eq.(5), which, as recognized already by Heisenberg, is 
more fundamental than the x-p commutator identity. 

2.4 Prom physics back to combinatorics 

To find a way out of the difficulties just described let us return back to the ex- 
pression < n \ O \ k > exp{iLLj{n,k)t). Suppose that the algebra of observables 
contains an identity element (operator). Then, by replacing O by this operator 
we obtain, < n \ k >=< n{t) \ k{t) > .This makes sense only if we require 
< n \ k >= constSnk- Clearly, we can always adsorb the constant into the 
definition of the scalar product. In this work, following [29], we suggest to re- 
place the basic commutators, Eq.(20), by the requirement of orthogonality. This 
requirement is compatible with the requirement that the operators describing 
observables are Hermitian whose eigenfunctions are mutually orthogonal. In- 
stead of operators whose explicit form is difiicult to obtain we shall focus our 
attention on the properties of orthogonal functions and, more generally, on the 
properties of orthogonal polynomials (e.g. s\, etc.). Development of theory of 
orthogonal polynomials of several variables in connection with quantum exactly 
solvable model systems is an active area of current research [38,46,47]. Such 
an approach makes sense since it is known [48] that all one- variable orthogonal 
functions of exactly solvable problems in quantum mechanics [49] are obtainable 
as various limiting cases of the Gauss-type hypergeometric functions^''. Follow- 
ing ideas by Aomoto, Orlik and Terrao demonstrate that the hypergeometric 
functions of multiple arguments (of which the Gauss-type is just a special case) 
are expressible in the form of period integrals^°. By the principle of comple- 
mentarity all many-body exactly solvable quantum mechanical problems should 
have hypergeometric functions of multiple arguments as eigenfunctions. The 
most important fact for our developments lies in the observation that when 

^^This fact will be discussed in detail in Section 3. 

^"Periods can be associated with the homology basis -different for different (algebraic) man- 
ifolds. Interested readers may consult either [48] or [50] for more details. 
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these functions become cigcnfunctions (as it is known in one component case), 
this produces orthogonal polynomials-different for different many- body quan- 
tum mechanical problems. This fact can be formulated as a problem : for a 
given set of orthogonal polynomials find the corresponding many-body oper- 
ator for which such a set of orthogonal polynomials forms a complete set of 
eigenfunctions. 

After these general remarks, we are ready to provide more concrete evi- 
dence that this is indeed the case. The symmetric group Sn has the following 
presentation in terms of generators Sj and Coxeter relations: 

si = 1, 

SiSj = SjSi for \i - j\ > 2, 
SiSi+iSi = Si+iSjSj+i. (21) 

If there is a set of n elements of whatever kind the generator Sj interchanges an 
element i with i -|- 1 so that si, ...,s„_i generate S'n. There are n! permutations 
in the set of n elements. If we assign the initial ordered state, then any other 
state can be reached by successful application of permutational generators to 
this state. The word w =SaiSa2 ■ --Sai (where the indices ai,...,ai represent a 
subset of the set of n — 1 elements) can be identified with such a state. Since 
one can reach this state in many ways, it makes sense to introduce the reduced 
word w whose length l(w) is minimal. We would like the generators of Sn to 
act on monomials = For this purpose, following Lascoux and 

Schiitzenberger [51] (L-S) we introduce an operator di via rule: 

:= . (22) 

Xi - Xi+i 

It acts on monomials such as in such a way that the generator Sj acting on 
the combination x'^'x'-^-^^ converts it into x'-'^^ x'^^-^. By construction, the action 
of this operator on monomial is zero if Ui = Oj+i, otherwise it diminishes the 
degree of the monomial by 1. In addition, the same authors introduce operators 

_^ ^ il^ilL^ (23) 

Xi - Xt+i 

and TTi = l-t-Tfj. Finally, being armed with such definitions, we can introduce an 
operator Di [p, q, r) = pdi + qiti + rsi'^^ , where p, q and r are some numbers. L-S 
demonstrate that such defined operator, while acting on monomials, obeys the 
braid- type relations (the 2nd and the 3-rd lines in Eq.(21)) while the relation 
sf = 1 is replaced by 

Df = qD^ + r{q + r). (24a) 

With constants p, q and r properly chosen, such a relationship defines the Hecke 
algebra Hn of the symmetric group 5„. Usually, it is written as 

= (1 - Q)Di + Q (24b) 



^^By doing so, the operators di,Tri and Sj become equivalent in the sense which we shall 
explain shortly. 
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with Q being some number (effectively playing the same role as q, r) . Hn 
should be considered as a deformation of iS'n.The rationale for its introduction 
lies in its direct connections with the knot and link theory so that quantum 
mechanics can be considered as some branch of this theory (Kholodenko 2006a). 
This fact will have its impact on quantization. To demonstrate this, following 
Kirillov Jr. [52], by relabeling earlier defined operator di as bij we reserve the 
notation di = g|- for the usual operator of differentiation. With its help we 
introduce the so called Dunkl operator T>i via 

Vi = di + kJ2bij (25) 

with k being some (known) constant. Such defined operator acts on monomials 

(polynomials). It possesses the property w'DiW~l='Dy,{i) Vi G Sn- Consider 
now the commutator \Pi,Vj\. Kirillov demonstrated that such a commutator is 
zero if bij satisfy the classical Yang-Baxter equation (GYBE) 

[bl2, fcis] + [&12, &23] + [&13, &23] = 0. (26) 

Alternatively, Eq.(26) can be taken as the definition for bij. This is facilitated by 
designing of the degenerate affine Hecke algebra (Cherednik 2005). The purpose 
of this algebra from the physical point of view is to introduce the Heiscnberg 
commutation rule Eq.(20) without reference to the B-S quanization prescription 
or to the (optical) sum rule described in Appendix A. Such an algebra is made 
up as a semidirect product of 5„ with the commutator algebra 

^i-^i — ^) XiSj — SjXi ^ j^j ~t~ -1) XiXj — XjXi^ (2Ta) 

where the constant h is playing essentially the same role as the Plank's constant 

h. From the above definitions it follows that Eq.(27a) is the discrete analog of 
the Heisenberg's commutation rule, Eq.(20). Furthermore, in view of the remark 
made after introduction of Di{p,q,r), it is possible to rewrite the commutator 
in Eq.(27a) in the equivalent form. This indeed was accomplished in the paper 
by Adin et al [53]. Hence, we can rewrite Eq.(27a) equivalently as 



XiQi OiXi-\-X — /i, OiXi - Xi-^idi — h , XiOj — OjXi Vi J, XiXj — XjXi^ 

(27b) 

where di should be understood in the sense of Eq.(22). At this point it is useful 

to introduce yet another operator Si = Si + hbi^i+i. It is designed in such a way 
that it obeys the braid relations: 

S1S2S1 = S2SiS2- (28) 

Furthermore, if now we define the operators R12 = sisi,i?23 = S2S2,^i3 = 
■S1-R23S1 = 52-^12*2, then the Eq.(28) becomes equivalent to the standard Yang- 
Baxter (Y-B) equation for Rij = 1 + hbij (or Rij ~ exp{hbij) for /i — > 0). 
Explicitly, we obtain: R12R13R23 = -^23-^13-^12- 
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All this discussion looks a bit formal at this point. Indeed, why to intro- 
duce the operator Pj? Why to be concerned about the commutator [Di,Dj]l 
What the Yang-Baxter equations have to do with the results of this paper? We 
would like to provide the answers to these questions now and in the following 
subsection. 

First, consider an equation 2?^/ = 0. It can be written alternatively as 




which is just the celebrated Knizhnik-Zamolodchikov (K-Z) equation^^. This 
means that: a) the operator Vi is effectively a covariant derivative (the Gauss- 
Manin connection in the formalism of fiber bundles) and, b) that the vanishing 
of commutator [DijVj] is just the zero curvature condition [54] essential for 
all known exactly intcgrablc systems. The question still remains: how 17^ in 
Eq.(29) is related to bij in Eq.(25)? The answer was found by Belavin and 
Drinfel'd [55]. In the simplest (rational) case we have bij{z) = as expected. 
More complicated trigonometric and elliptic cases found by Belavin and Drin- 
fel'd are summarized in the book by Etingof with collaborators [56]. From the 
references just provided, it should be clear that since solutions of the K-Z equa- 
tions are expressible; in terms of hypergeometric functions of single and multiple 
arguments, all examples of exactly solvable quantum mechanical problems (in- 
cluding those involving the Dirac equation, and, hence, the spin) found in the 
textbooks on quantum mechanics are covered by the formalism we have just 
described. In the next subsection we would like to illustrate these results by 
concrete physical examples taken from current physical literature. 



2.5 Latest developments in atomic physics illustrating gen- 
eral principles 



In the review paper by Tanner et al [57] as well as in the book by Cvitanovic'[58] 
it is explained in detail that in order to calculate the He spectrum it is sufficient 
either: a) to consider the classical dynamics of two electrons and the nucleus on 
the line and to use this information in the semiclassical trace formula producing 
very accurate results for the spectrum or, b) to restrict quantum mechanical cal- 
culations to the spherical approximation (the so called s- wave approximation) 
in order to arrive at the exactly solvable radial Schrodingcr-typc equation for 
two electrons and massive nucleus[59,60] producing very reasonable results for 
the He spectrum. To these achievements we would like to add those by Svidzin- 
sky et al [10,11] and Muravski and Swidzinsky [12] where the same type or even 
better accuracy for He and other atoms and diatomic molecules is obtained 
using the so called d-scaling. In this method the multielectron Schrodinger 
equation is analyzed in various dimensions. Upon proper rescaling, the limiting 

^^In fact, in general case [48] scalar function /(z) is replaced by the vector function f (z).This 
fact should be kept in mind in actual calculations. 
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case: oo, is reduced again to the exactly solvable radial- type multielectron 
equation which in the present case becomes classical equation considered al- 
ready by Bohr in 1913. Thus again, the zeroth order exactly solvable problem 
is one dimensional. Corrections in powers of 1/d are easily calculable produc- 
ing results which compare extremely well with much more cumbersome (and 
time consuming) Hartree-Fock type calculations. To this list of examples it is 
appropriate to add work by Ostrovsky and Prudov [61] which uses essentially 
the same averaging and perturbation methods as developed in celestial mechan- 
ics [62] superimposed with the Bohr-Sommerfeld quantization prescription. All 
examples discussed in this subsection wore done without theoretical guidance 
(other than the proof of the existence of minimizers for Bohr-type functionals 
[13]. The theoretical framework developed in this section naturally explains 
why these results are actually working so well. This framework sets up the 
stage for developing applications of these results to celestial mechanics to be 
discussed in the rest of this paper. 

3 Space, time and space- time in classical and 
quantum mechanics 

3.1 General comments 

If one contemplates quantization of dynamics of celestial objects using tradi- 
tional textbook prescriptions, one will run into myriad of small and large prob- 
lems immediately. Unlike atomic systems in which all electrons repel each other, 
have the same masses and are indistinguishable, in the case of, say, Solar Sys- 
tem all planets (and satellites) attract each other, have different masses and 
visibly distinguishable. Besides, in the case of atomic systems the Planck con- 
stant h plays prominent role while no such a role can be given to the Planck 
constant in the sky. The only thing which remains in common between both 
atomic and celestial dynamic systems is the existence of stable closed orbits. In 
the previous section we demonstrated that this fact is absolutely essential for 
quantization. Nevertheless, the formalism developed thus far resembles more 
the existence theorem in mathematics rather then the actual manual describing 
the computational protocol. The task now lies in developing necessary construc- 
tive steps leading to actual implementation of general principles. This task is 
accomplished below and in the following section 

3.2 Space and time in clcissical and quantum mechanics 

Although celestial mechanics based on the Newton's law of gravity is consid- 
ered to be classical (i.c.nonquantum), with such an assumption one easily runs 
into serious problem. Indeed, such an assumption implies that the speed with 
which the interaction propagates is infinite and that the time is the same every- 
where. Wether this is true or false can be decided only experimentally. Since 
at scales of our Solar System one has to use radio signals to check correctness 
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of Newton's celestial mechanics, one is faced immediately with all kind of wave 
mechanics effects such as retardation, the Doppler effect, etc. Because of this, 
the measurements are necessarily having some error margins. The error margins 
naturally will be larger for more distant objects. Accordingly, even at the level 
of classical mechanics applied to the motion of celestial bodies we have to deal 
with certain inaccuracies similar in nature to those in atomic mechanics. To 
make formalisms of both atomic and celestial mechanics look the same one has 
to think about the space, time and space-time transformations already at the 
level of classical mechanics. 

We begin with observation that in the traditional precursor of quantum 
mechanics-the Hamiltonian mechanics-the Hamiltonian equations by design re- 
main invariant with respect to the canonical transformations [63]. That is if 
sets {qi} and {pi} represent the "old" canonical coordinates and momenta while 
Qi = Qi{{qi}, {Pi}) and Pi = Pi{{qi}, {ft}), i = 1- N, represent the "new" set 
of canonical coordinates and momenta, the Hamiltonian equations in the old 
variables given by 

dH ^ . dH 
Qi = and pi = -— (30) 
dpi dqi 

and those rewritten in " new" variables will have the same form. Here we used 
the commonly accepted notations, e.g. Qi = -^qi , etc. Quantum mechanics uses 
this form-invariance essentially since the Poisson brackets introduced in Eq.(9) 
by design will also have the same form in terms of both "old" and "new" 
canonical variables. 

We would like to complicate matters by investigating the possibility of the 
"canonical " time changes in classical mechanics. Fortunately, such a possibility 
was explored to a some extent already. This is described in the monograph by 
Pars [63] thus making our task considerably simpler. For the sake of space, we 
refer our readers to pages 535-540 of this monograph. Furthermore, following 
Dirac [64], we notice that a good quantization procedure should always begin 
with the Lagrangian formulation of mechanics since it is not always possible to 
make a transition from the Lagrangian to Hamiltonian form of mechanics (and, 
thus, to quantum mechanics) due to the presence of some essential constraints 
( typical for mechanics of gauge fields, etc.). Hence, we also begin with the 
Lagrangian functional C = C{{qi}, {qi})- The Lagrangian equations of motion 
can be written in the form of Newton's equations given by Pi = Fi, where the 
generalized momenta pi are given by Pi = 5C/5qi and the generalized forces Fi 
are given by Fi = —5C/5qi. In the case if the total energy E is conserved, it is 
possible instead of "real" time t to introduce the fictitious time 9 via relation 
dt = u{{qi})d6 where the function u{{qi]) is assumed to be nonnegative and is 
sufficiently differentiable with respect to its arguments. At this point we can 
enquire if Newton's equations can be written in terms of new time variable so 
that they remain form- invariant. To do so, following Pars, we must: a) to 
replace C by uC, b) to replace qi by q[ /u, where q'i=2eqii c) to rewrite the 
new Lagrangian in terms of such defined new time variables and, finally, d) to 
obtain Newton's equations according to the described rules, provided that now 
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wc have to use p[ instead of pi. In the case if the total energy of the system is 
conserved, we shall obtain back the same form of Newton's equations rewritten 
in terms of new variables. This means that by going from the Lagrangian 
to Hamiltonian formalism of classical mechanics we can write the Hamilton's 
Eq.(30) in which the dotted variables are replaced by primed. These arguments 
demonstrate connections between space and time already at the level of classical 
mechanics. Situation here is similar to that encountered in thermodynamics 
where instead of absolute temperature one can use any nonegative function 
of absolute temperature as new temperature. Using these arguments we notice 
that since the temperature is conjugate to energy in thermodynamics, the time is 
conjugate to energy in mechanics and, accordingly, in quantum mechanics. This 
means that for the nondissipative (i.e. energy conserving) Hamiltonian system^^ 
the Hamiltonian equations of motion, Eq.(30), will remain form- invariant if we 
replace the Hamiltonian H by some nonnegative function f{H) while changing 
time t to time 6 according to the rule d6/dt = df{H)/dH \h=e- Such a change 
will affect the quantum mechanics where now the Schrodinger's equation 

ih^^ = (31a) 
at 

is to be replaced by 

zfi^* = /(iJ)vl>. (31b) 

With such an information at our hands, we would like to discuss the extent to 
which symmetries of our (empty) space-time affect dynamics of particles " living" 
in it. 

3.3 Space-time in quantum mechanics 
3.3.1 General comments 



Use of group-theoretic methods in quantum mechanics had began almost imme- 
diately after its birth. It was initiated by Pauli inl926. He obtained a complete 
quantum mechanical solution for the Hydrogen atom employing symmetry argu- 
ments. His efforts were not left without appreciation. Our readers can find many 
historically important references in two comprehensive review papers by Bander 
and Itzykson [65]. In this subsection we pose and solve the following problem: 
Provided that the symmetry of (classical or quantum) system is known, will 
this information be sufficient for determination of this system uniquely? Below, 
we shall provide simple and concrete examples illustrating meaning of the word 

^^It should be kept in mind that the concept of nondissipativity is actually of quantum 
origin (e.g. recall superconductors or superfiuids). In classical mechanics such a concept is 
just a convenient idealization similar to the notion of a material point in Newton's mechanics 
or the notion of thermodynamics when it is applied to the real heat engines, etc. The truly 
nondissipative mechanical systems thus should behave quantum mechanically. This observa- 
tion provides the hint that some stable motions in our Solar System are of quantum nature. 
In view of Eq. (6) this option maJses sense. 
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" determination" . In the case of quantum mechanics this problem is known as 
the problem about hearing of the " shape of the drum" . It was formulated by 
Mark Kac [66]. The problem can be formulated as follows. Suppose that the 
spectrum of the drum is known, will such an information determine the shape 
of the drum uniquely? The answer is " No" [67] . Our readers may argue at this 
point that nonuniqueness could come as result of our incomplete knowledge of 
symmetry or, may be, as result of the actual lack of true symmetry (e.g. the 
Jahn- Teller effect in molecules, etc. in the case of quantum mechanics). These 
factors do play some role but they cannot be considered as decisive as the basic 
example below demonstrates. 

3.3.2 Difficulties with the correspondence principle for Hydrogen 
atom 

In this subsection we even do not use arguments by Kac. Since our arguments 
are straightforward, they are more intuitively appealing. We choose the most 
studied case of Hydrogen atom as an example. 

As it is well known, the classical mechanical problem about motion of the 
particle in centrally symmetric field is planar and is exactly solvable for both 
the scattering and bound states [63,68]. The result of such a solution depends 
on two parameters: the energy and the angular momentum. The correspon- 
dence principle formulated by Bohr is expected to provide the bridge between 
the classical and quantum realities by requiring that in the limit of large quan- 
tum numbers the results of quantum and classical calculations for observables 
should coincide. Appendix A provides a good example of such kind of thinking. 
However, this requirement may or may not be possible to implement. It is vio- 
lated already for the Hydrogen atom. Indeed, according to the naive canonical 
quantization prescriptions, one should begin with the classical Hamiltonian in 
which one has to replace the momenta and coordinates by their operator analogs. 
Next, one uses such constructed "quantum" Hamiltonian in the Schrodinger's 
equation, etc. Such a procedure breaks down at once for the Hamiltonian of 
Hydrogen atom since the intrinsic planarity of the classical Kepler's problem 
is entirely ignored thus leaving the projection of the angular momentum with- 
out its classical analog. Accordingly, the scattering states of Hydrogen atom 
and the classical mechanically obtained Rutherford's formula obtained for pla- 
nar configurations are reproduced quantum mechanically (within the 1st Born 
approximation) using the 3-d Schrodinger's equation ! Thus, even for the Hydro- 
gen atom the classical and the quantum (or, better, pre quantum) Hamiltonians 
do not match thus formally violating the correspondence principle. Evidently, 
semiclassically we can only think of energy and the angular momentum thus 
leaving the angular momentum projection unobserved. Such a "sacrifice" is 
justified by the agreement between the observed and predicted Hydrogen atom 
spectra and by use of Hydrogen-like atomic orbitals for multielectron atoms. 
Although, to our knowledge, such a mismatch is not mentioned in any of the 
students textbooks on quantum mechanics, its existence is essential if we are in- 
terested in applications of quantum mechanical ideas to Solar System dynamics. 
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In view of such an interest, we would like to reconsider traditional treatments 
of Hydrogen atom, this time being guided only by symmetry considerations. 
This is accomplished in the next subsection. 

3.3.3 Emergence of the SO(2,l) symmetry group 

In April of 1940 Jauch and Hill [69] published a paper in which they studied the 
planar Kepler problem quantum mechanically. Their work was stimulated by 
earlier works by Fock of 1935 and by Bargmann of 1936 in which it was shown 
that the spectrum of bound states for the Hydrogen atom can be obtained by 
using representation theory of SO (4) group of rigid rotations of 4-dimensional 
Euclidean space while the spectrum of scattering states can be obtained by using 
the Lorentzian group S0(3,l). By adopting results of Fock and Bargmann to 
the planar configuration Jauch and Hill obtained the anticipated result. In 
the planar case one should use S0(3) group for the bound states and S0(2,l) 
group for the scattering states. Although this result will be reconsidered almost 
entirely, we mention about it now having several purposes in mind. 

First, we would like to reverse arguments leading to the final results of Jauch 
and Hill in order to return to the problem posed at the beginning of this section. 
That is, the fact that the Kepler problem is planar (due to central symmetry 
of the force field) and the fact that the motion is restricted to the plane and 
takes place in (locally) Lorentzian space-time are the most general symmetry 
constraints imaginable. Thus, the fact that the Lorentz S0(2,l) group is related 
to the spcctriim of Kepler problem should be anticipated. Nevertheless, the 
question remains: is Kepler's problem the only one exactly solvable classical 
and quantum mechanical problem associated with the S0(2,l) group? Below 
wc demonstrate that, unfortunately, this is not the case. In anticipation of such 
negative result, we would like to develop our intuition by using some known 
results from quantum mechanics. 

3.3.4 Classical-quantum correspondence allowed by SO(2,l) symme- 
try: a gentle introduction 

For the sake of space, we consider here only the most generic (for this work) 
example in some detail: the radial Schrodinger equation for the planar Kepler 
problem with the Coulombic potential. It is given by^^ 

2/x dp'^ pap p'^ p 

Here |m| = 0, 1, 2, .. is the angular momentum quantum number as required. For 
S < it is convenient to introduce the dimensionless variable x via p = ax and 
to introduce the new wave function: ■^(p) = y/p^'(p). Next, by the appropriate 

^*The rationale for discussing the Coulombic potential instead of gravitational will be fully 
explained in the next section. 
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choice of constant a and by redefining ij){p) as ^(p) = 7a;2+l'"l exp(— 
where y = 7X, -7^ = '^^a'^,a = the following hypergeometric equation 

can be eventually obtained: 

{y^. + 2[|m| + ^ - + 2[1 - \m\ - '-]] ^(y) = 0. (33) 

Formal solution of such an equation can be written as ip{y) = !F{—A{m), B{m), y), 
where !F is the confluent hypergeometric function. Physical requirements im- 
posed on this function reduce it to a polynomial leading to the spectrum of 
planar Kepler problem. Furthermore, by looking into standard textbooks on 
quantum mechanics, one can easily find that exactly the same type of hyperge- 
ometric equation is obtained for problems such as one-dimensional Schrodinger's 
equation with the Morse-type potential, three dimensional radial Schrodinger 
equation for the harmonic oscillator^^ and even three dimensional radial equa- 
tion for the Hydrogen atom^''. Since the two-dimensional Kepler problem is 
solvable with help of the representations of S0(2,l) Lorentz group, the same 
should be true for all quantum problems just listed. That this is the case is 
demonstrated, for example, in the book by Wybourne [70]. A sketch of the 
proof is provided in Appendix B. This proof indicates that, actually, the dis- 
crete spectrum of all problems just listed is obtainable with help of S0(2,l) 
group. The question remains: if the method outlined in Appendix B provides 
the spectra of several quantum mechanical problems listed above, can we be 
sure that these arc the only exactly solvable quantum mechanical problems as- 
sociated with the S0(2,l) Lorentz group? Unfortunately, the answer is "No". 
More details are given below. 

3.3.5 Common properties of quantum mechanical problems related 
to SO(2,l) Lorentz group 

In Appendix B we provide a sketch of the so called spectrum-generating algebras 
(SGA) method producing the exactly solvable one-variable quantum mechani- 
cal problems. In this subsection we would like to put these results in a broader 
perspective. In particular, in Section 2 we demonstrated that all exactly solv- 
able quantum mechanical problem should involve hypergeometric functions of 
single or multiple arguments. We argued that the difference between different 
problems can be understood topologically in view of the discussed relationship 
with braid groups. On another hand, obtained results, even though rigorous, are 
not well adapted for immediate practical use. In this regard more useful would 
be to solve the following problem: For a given set of orthogonal polynomials 

^^That is, V{x) = A{exp{-2ax) — 2exp{—ax)). 
2*'That is, V{r) = ^ Br'^. 

27That is, y(r) = 4 - -■ 
r 
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find the corresponding many-body operator for wliicli sucli a set of ortliogonal 
polynomials forms a complete set of eigenfunctions. At the level of orthogonal 
polynomials of one variable relevant for all exactly solvable two-body problems 
of quantum mechanics, one can think about related problem of finding all poten- 
tials in the one-dimensional radial Schrodinger equation, e.g. Eq.(B.l), leading 
to the hypergeometric-type solutions. Such a task was accomplished by Natan- 
zon [71]. Subsequently, his results were reinvestigated by many authors with 
help of different methods, including SGA. To our knowledge, the most complete 
recent summary of the results, including potentials and spectra can be found in 
the paper by Levai [72]. Even this (very comprehensive) paper does not cover 
all aspects of the problem. For instance, it does not mention the fact that these 
results had been extended to relativistic equations such as Dirac and Klein- 
Gordon for which similar analysis was made by Cordero with collaborators [73] 
. In all cited cases (relativistic and non relativistic) the underlying symmetry 
group was S0(2,l). The results of Appendix B as well as all other listed ref- 
erences can be traced back to the classically written papers by Bargmann [74] 
and Barut and Fronsdal [75] on representations of S0(2,l) Lorentz group. Fur- 
thermore, the discovered connection of this problematic with supersymmetric 
quantum mechanics [76,77] can be traced back to the 19th century works by 
Gaston Darboux [72]. 

Summarizing, established in Section 2 rigorous connections between exactly 
solvable two-body quantum mechanical problems and hypergeomctric functions 
and, by complementarity principle, between exactly solvable many body prob- 
lems and hypergeomctric functions of many arguments are consequences of the 
locally Lorentzian group structure of our space-time. Such a structure allows 
many but not infinitely many exactly solvable problems to exist. The fact that 
planar S0(2,l) is sufficient to cover all known exactly solvable two-body cases 
(instead of the full S0(3,l) Lorentz group!) is quite remarkable. It is sufficient 
for the purposes of this work but leaves open the question : Will use of the full 
Lorentz group lead to the exactly solvable quantum mechanical problems not 
accounted by S0(2,l) group symmetry? This topic will be discussed in Section 
5. In the meantime, we would like to address the problem of quantization of 
Solar System dynamics using results of Sections 2 and 3. This is done in the 
next section 

4 Quantum celestial mechanics of Solar System 
4.1 General remarks 

We begin this section by returning back to Eq.(6) once again. Based on previous 
discussions, this equation provides us with the opportunity to think seriously 
about quantum nature of our Solar System dynamics. Nevertheless, such an 
equation reveals only one aspect of quantization problem and, as such, provides 
only a sufficient condition for quantization. The necessary condition in atomic 
and celestial mechanics lies in the nondissipativity of the dynamical systems 
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in both cascs^^. Recall that Bohr introduced his quantization prescription to 
avoid dissipation caused by the emission of radiation by electrons in orbits in 
general position. As we demonstrated previously, new quantum mechanics have 
not explained absence of dissipation for stationary Bohr's orbits^^. In fact, as 
our analysis of Heisenberg's work(s) indicates, new quantum mechanics have 
not added a single new element to the old atomic mechanics in terms of the new 
issues to be considered. 

In the nutshell, new quantum mechanics provided a convenient computa- 
tional scheme for dealing with otherwise purely mechanical problems involv- 
ing accidental degeneracy (that is resonances). By doing so, it made no at- 
tempt at explaining (using known results from mechanics and electrodynam- 
ics) the nondissipativity. Nevertheless, the phenomenon of nondissipativity was 
explained quite convincingly in the case of superconductivity and superfluidity 
later on. Thanks to these intrinsically quantum phenomena, we can be sure that 
quantum mechanics did capture some truth. Regrettably, only some since, as 
we discussed in Section 3.2.2, even for the most studied case of Hydrogen atom 
the task of establishing the correspondence between the classical and quan- 
tum models of Hydrogen atom is nontrivial. The symmetry (and supersymme- 
try) arguments of Section 3 based on locally Lorentzian space-time structure as 
well as the combinatorial arguments of Section 2 simplified task of establishing 
the quantum-classical correspondence considerably. This happened because of 
firmly established finite number of exactly solvable quantum mechanical prob- 
lems allowed by the Lorentzian-type symmetry whose spectra are known and 
documented. These facts allow us to think seriously about quantization of Solar 
System dynamics. 

4.2 From Laplace to Poincare' and Einstein 

Before discussing this issue in some detail, we still need to make several remarks. 
First, although superficially classical Hamiltonians for Coulombic and Newto- 
nian potentials look almost the same, the naive textbook-style quantization will 
immediately run into major problems. For one thing, all electron masses are the 
same while all planetary/satellite masses are different. For other thing, filling 
of atoms by electrons is controlled by the electric charge of the nucleus so that 
stable atoms/molecules are electrically neutral. Apparently, no such restriction 
exists for the system of gravitating bodies. Next, apparent violation of pla- 
narity of Hydrogen atom treated at the level of classical mechanics is justified 
by the fact that the angular momentum projection does play an important role 
in chemistry. As far as we can see, nothing of that sort exists in the sky. 

To deal with the mass differences for planetary systems we have to recall 
some facts from general relativity. We shall restrict ourself only by some illus- 
trative examples meant to provide some feeling of problems we would like to 
discuss. To this purpose we would like to make some comments on the classical 

^^E.g. see the paper by Goldreich [25] mentioned in Section 1.2. 

^^At the level of old Bohr theory absence of dissipation at the stationary Bohr orbit was 
explained by Boyer [78]. Subsequently his result was refined by Puthoff [79]. 
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mechanical treatment of Kepler problem in representative physics textbooks, 
e.g. read [68, 80]. Such treatments tend to ignore the equivalence principle essen- 
tial for the gravitational Kepler problem and nonexistent for the Coulomb-type 

problems. This causes some significant inaccuracies to emerge. Specifically, ac- 
cording to Vol.2 of famous Landau-Lifshitz course in theoretical physics [81] if 

2 

we take L= — ^ nup as the Lagrangian for a particle in gravitational field 

(represented by a local potential (p),the Lagrangian (Newtonian) equations of 
motion can be written as 

V = -V(p (34) 

so that the mass drops out of this equation making it possible to think about 
such an equation as an equation for geodesic in pseudo-Riemannian space. This 
observation had lead Einstein to full development of general relativity theory. 
By noticing that Newton's equation makes sense only for material points ( that 
is for idealized formally nonexisting objects), the same must be true for Eq.(34). 
Hence, as such it is valid only for the well localized point-like objects. Using 
such idealized model, we need to discuss briefly the 2-body Kepler problem for 
particles with masses mi and m2 interacting gravitationally. The Lagrangian 
for this problem is given by 

mi .2 , m2 .2 , mim2 , , 

C = —r^ + —r2+^- — . (35a) 

2 2 |ri-r2| 

Introducing the center of mass and relative coordinates via miri -|- m2r2 = 
and r = ri — r2 , the above Lagrangian can be rewritten as 

/i.2 , mim2 mim2 P (mi-hm2), 
i'=— r +7 — — ; — = + 7 r-, (35b) 

2 ' |r] mi + m2^ 2 ' |r| " ^ ' 

where, as usual, we set u = "^^^ .The constant can be dropped and, 

' ' ^ mi-t-m2 mi+m2 ' 

after that, instead of the geodesic (34) we obtain the equation for a fictitious 
point-like object of unit mass moving in the field of gravity produced by the 
point-like body of mass mi H- m2. Clearly, in general, one cannot talk about 
geodesies in this case. Nevertheless, as it is usually done, if, say, mi » m2 (as 
for the electron in Hydrogen atom or for the Mercury rotating around Sun) one 
can with very good accuracy discard mass m2 thus obtaining an equation for 
a geodesic. Such an approximation was indeed made by Einstein in his major 
work on general relativity [82] in which he ignored the mass of Mercury entirely 
when making his calculations of the perihelium shift for this planet. More recent 
results [83] show that such an approximation is expected to be quite satisfactory 
for other planets of our Solar System^°. With the exception of Pluto-Charon 
system, where /U2 = m2/(mi-|-m2) is of order 10~^, and the Earth-Moon system, 
where /X2 is of order 10~^, all other planet-satellite and Sun- planet pairs have 
/X2 of order 10~^ and less [26] so that use of geodesies is justifiable physically. 
Mathematically, however, this is not quite the case yet since, even in the case 



■'E.g. read Box 40.3 of this reference as well as pages 1126-1129. 
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of Mercury considered by Einstein, it is necessary to prove that influence of the 
rest of planets of Solar System on its motion can be ignored (as well as finite 
size of the Sun, etc.). 

The task of proving that motion of planets can be well approximated by 
geodesies can be traced back to works by Laplace on celestial mechanics. Lec- 
ture notes by Moser and Zehnder [84] contain accessible discussion of Laplace's 
works^"*^ to which we refer our readers for details. In short, Laplace was interested 
in dynamics of the planar 4-body problem using Jupiter and its 3 satellites: 
lo, Europe and Ganymede, as an example. He noticed that the motion of these 
satellites obeys the resonance condition and he was able to reproduce this mo- 
tion analytically by ignoring satellite masses (just like in Eq.(35b), but beginning 
with the full 4-body problem initially !). Under these conditions, gravitational 
interactions between satellites can be neglected so that the motions become 
completely decupled but subject to the resonance condition. Furthermore, to 
study stability of such resonance motions Laplace (and Lagrange) assumed that 
the actual (Lagrangian) motions^^ of satellites oscillate about the respective 
stable orbits of these satellites. Thus, effectively, Laplace and Lagrange were 
considering the effects of general relativity and quantum mechanics long before 
these disciplines have been officially inaugurated. In their lectures, Moser and 
Zehnder also provide references to works by Poincare' and de Sitter on further 
refinements of Laplace's results. Although according to Arnol'd et al [18] the 
extension of work by Laplace to the full n+1 body planar problem was given 
in the monograph by Charlier [85] , rigorous mathematical proofs have been ob- 
tained only quite recently by Fejoz [86] and Biasco et al [87]. To realize the 
difficulties in providing such a proof it is sufficient, following Poincare' [88], to 
demonstrate that the results of massless limit considered by Laplace will re- 
main practically unchanged if the satellites would have some small but finite 
masses (so that they interact with each other !). Such a philosophy lies at the 
heart of KAM theory used and improved in the works by Fejoz [86] and Biasco 
et al [87]. 

Even with these proofs available, one should take into account that, in view 
of experimental limitations, Newton's law of gravity should amended rigor- 
ously speaking. This is so taking into account the finite speed of propagation of 
gravitational interaction as well as the fact that all observations are made with 
some kind of light/radio sources causing retardation, Doppler and other effects. 
Thus, taking into account experimental conditions, the traditional classical 
mechanics description of celestial motions becomes replaced by that encoun- 
tered in quantum mechanics where one has to use probabilities to account for 
incompleteness of available information. Furthermore, the above proofs do not 
account for dissipation effects playing major stabilizing role in both atomic and 
quantum celestial mechanics. 

If we assume that the motion of bodies indeed takes place on geodesies then, 
formally, there are no interactions and the local time becomes proper time. In 

31 E.g. read pages 102-120. 

32E.g. see Arnol'd et al [18], page 261. 
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the case of, say, binary stars of comparable masses one cannot use geodesies for 
description of their relative motion'^'^ so that one is confronted with the problem 
of matching the Einsteinian gravity with its Newtonian limit as discussed in the 
paper by Einstein, Infeld and Hoffmann [89] In all other theories of gravity, 
including the Brans-Dicke- Jordan's theory, there are substantial departures from 
the geodesic motion. Details can be found on pages 1127-1129 of the book by 
Misner et al [83]. 

Clearly, the difficulties of explaining motions using classical mechanics of 
n+1 body problem are such that the assumption about truly geodesic motion 
looks suspicious. But these results are based on Newtonian mechanics which 
by design do not account for dissipation and retardation effects. The facts 
we just mentioned also complicate the choices between different (alternative) 
theories of gravity. Hence, it is clear that at the present state of our knowledge 
the ultimate choice between competing theories can be made only based on 
additional information. Such an information is supplied, in part, in this work 
where uses of historical analogies between the quantum (atomic) and celestial 
mechanics provide some helpful guidance. For this purpose we compiled our 
Table 3 prior to actual computations. 



Table 3 



\Type of mechanics 
Properties 


Quantum atomic 

mechanics 


Quantum 

celestial mechanics 


Dissipation (type of)\ 
(yes\no)\on stable orbits 


electromagnetic 
friction\no\ 
Bohr orbits 


tidal friction 
\no\Einstein's geodesies 


Accidental degeneracy\ 
(yes\no)\origin 


yes\Bohr-Sommerfeld condition 


yes\closure of the Lagrangian orbits 


Charge neutrality 


yes 


no (but see below) 


Masses 


electrons having 
the same masses 


(up to validity of the 
equivalence principle) 
masses are the same 


Minimal symmetry group 


S0(2,l) 


S0(2,l) 


Correspondence 
principle 


occasionally violated 


occasionally violated 


Discrete spectrum: 
finite or infinite\reason\ 
Pauli principle(yes\no) 


finite and infinite\ 
charge neutrality\ 
yes 


finite\ 

no charge neutrality\ 
yes 



•^■'This case was discussed in papers by Einstein, Infeld and Hoffmann [89] and Robertson 
[90] with the outcome that it is possible to describe gravitational field outside such a binary 
system in terms of geodesies. This leaves open the question of dynamical stability of such 
binaries since their motion is controlled by the Newton's equations of motion. In view of 
the effects of tidal friction, which should be quite appreciable in this case, the dynamics of 
such binaries should be most likely unstable. For such systems one can safely neglect friction 
caused by the emission of gravitational waves since these are effects of fifth order in c~^(e.g. 
read Landau and Lifshitz [81], paragraph 106). 

^*See also Section 6 below. 
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Details related to this table are discussed further below. 

4.3 Celestial spectroscopy, the Titius-Bode law of plane- 
taiy distances and quantum celestial mechanics 

The atomic spectroscopy was inaugurated by Newton, in the second half of 
17th century. As wc discussed in the Introduction, the results of atomic and 
molecular spectroscopy were used by Bohr in essential way resulting in the birth 
of quantum mechanics. The celestial spectroscopy was inaugurated by Titius 
in the second half of 18th century and become famous after it was advertised 
by Johann Bode, the Editor of the " Berlin Astronomical Year-book" in the late 
18th century. The book by Nieto [91] provides extensive bibliography related 
to uses and interpretations of the Titius-Bode (T-B) law up to second half of 
20th century. Unlike the atomic spectroscopy, where the observed atomic and 
molecular spectra were expressed using simple empirical formulas which were 
(to our knowledge) never elevated to the status of "law", in celestial mechanics 
the empirical T-B formula 

r„ = 0.4-h03. • 2", n = -00,0,1,2, 3,... (36) 

for the orbital radii (semimajor axes) of planets acquired the status of a law in 
the following sense. In the case of atomic spectroscopy the empirical formulas 
used for description of the atomic/molecidar spectra have not been used (to 
our knowledge) for making predictions. Their purpose was just to describe in 
mathematical terms what had been observed. Since the T-B empirical formula 
for planetary distances was used as the law, it was used in search for planets not 
yet discovered. In such a way Ceres, Uranus, Neptune and Pluto were found 
[92]. However, the discrepancies for Neptune and Pluto were much larger than 
the error margins allowed by the T-B law^^.This fact divided the astronomical 
community into "believers" and "atheists" regarding to the meaning and uses 
of this law. Without going into historical details, we would like to jump to the 
very end of the Titius-Bode story in order to use its latest version which we 
found in the paper by Neslusan [93] who, in turn, was motivated by the work of 
Lynch [94]. Instead of Eq.(36) these authors use another empirical power law 
dependence 

r„ = roB", n = l,2,3,..,9. (37) 

For planets (except Pluto and including the asteroid belt) Neslusan obtained^^ 
ro(au) = 0.203 and B = 1.773 with the rms deviation accuracy of 0.0534^'^. 
Analogous power law dependencies were obtained previously in the work by 
Dermott [95] for both planets and satellites of heavy planets such as Jupiter, 
Saturn and Uranus. 

•^^ Chapter 10 of [92] provides a very lively account of the present knowledge about various 
objects "living" in the Solar System. 
^®In astronomical units (to be defined below). 

^'^This result gives for the Earth in astronomical (au) units the result rs ~ 1.13. Much 
better result is obtained in case if we choose B = 1.7. In this case we obtain: rs ~ .997339. 
Lynch (2003) provides B = 1.706 and ro = 0.2139. 
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It should be noted that because of noticed discrepancies the attempts were 
made to prove or disprove the Titius-Bode law by using statistical analysis, e.g. 
see papers by Lynch [93] and Hayes and Tremaine [96], with purpose of finding 
out to which extent the observed dependencies can be considered as non ac- 
cidental. In view of Heisenberg quantization (honeycomb) condition, Eq.(5.a), 
it should be obvious by now that whatever distribution of frequencies can be 
measured, it can, in principle, lead to quantization. In principle, because to 
implement this in practice requires to identify possible models and the Hamil- 
tonians for these models as we discussed extensively in the previous sections. 
Hence, in the present case we arc confronted with exactly the same task. To 
move forward some historical analogies are helpful at this time. 

When Bohr was analyzing the data for He atom (Table 1) he had in mind 
a model of He made of two independent electrons rotating around the same 
nucleus. As results of Section 1 indicate, such an approximation produced 
quite reasonable results. Clearly, when dealing with dynamics of Solar System, 
one would like to follow the same philosophy. That is to assume first that 
the planets are noninteracting and move along the geodesies independently. In 
the case of atomic mechanics it was clear from the beginning that such an 
approximation should sooner or later fail even though it works well in some 
cases. For exactly the same reasons it is rather naive to expect that the T-B 
law makes always sense. Rather, it makes sense to assume that it works for as 
long as the assumption of noninteracting planets moving on geodesies can be 
checked quantum mechanically. Furthermore, the nonexisting electroneutrality 
in the sky provides strong hint that the T-B law must be of very limited use 
since the number of discrete levels for gravitating systems should be always 
finite. Otherwise we would observe the countable infinity of satellites around 
Sun or heavy planets which is both observationally and physically wrong. In 
the literature one can find many attempts at quantization of Solar System using 
sttandard prescriptions of quantum mechanics'^*. Because of this, restrictions 
on the number of allowed discrete levels cannot be made. 

In the present case, to facilitate matters, we would like to make several 
additional observations. First, we have to find the analog of the Planck constant. 
Second, we have to have some mechanical model in mind to make our search 
for correct answer meaningful. To accomplish the first task, we have to take 
into account the 3-rd Kepler's law. In accord with Eq.(35b), it can be written 
47r2 

as r^/T^ — r- In view of arguments presented in previous subsection, 

7(-M + to) 

we can safely approximate this result by 4:Tt^/jM, where M is the mass of 
Sun. For the purposes of this work, it is convenient to restate this law as 
3Znr„ — 21nT„ = In47r^/7M = const Below, wc choose the astronomical 
system of units in which Att'^ /jM = 1. By definition, in this system of units we 
have for the Earth: rs = T3 = 1. 

Consider now the Bohr result, Eq.(4), and take into account that E = 

*® Since this work is not a review, we do not provide references to papers whose results do 
not affect ours. 
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hui = ^.Therefore, Bohr's result can be conveniently restated as w{n,m) = 

u){n) — ti)(m). Taking into account Eq.s(4),(31b),(37) and the third Kepler's law 
we obtain: ^ 

w(n,m) = ^(ncln^ — mc In A), (38) 

cAnA 

where the role of Planck's constant is played now by clnA where A = Bi and 
c is some constant which will be determined selfconsistently below^^. 

At first, one may think that what wc obtained is just a simple harmonic 
oscillator spectrum. After all, this should come as not too big a surprise since in 
terms of the action -angle variables all exactly integrable systems are reducible 
to the sets of harmonic oscillators. This result is also compatible with the 
results of Appendix B. The harmonic oscillator option is physically undesirable 
in the present case since for gravitating systems the charge neutrality constraint 
cannot be imposed, e.g. sec Table 3. Evidently, allowing such a spectrum is 
equivalent to the correctness of the T-B law. But it is well known that this law 
is not working well for larger numbers. In fact, it would be extremely strange 
should it be working in this rc\Q,ime since the; total mass of all harmonically 
bound planets could potentially become infinite. 

To make a progress, we have to use the 3rd Kepler's law once again, i.e. we 
have to take into account that in chosen astronomical system of units 3Znr„ = 
2 lnr„. In view of the arguments just presented, a quick look at Eq.s B(13),(14) 
suggests that the underlying mechanical system is likely to be associated with 
that for the Morse potential. The low lying states of such a system cannot be 
distinguished from those for the harmonic oscillator. However, this system does 
have only a finite number of energy levels which makes sense physically. The 
task remains to connect this system with the planar Kepler's problem. Although 
in view of results of Appendix B such a connection does indeed exist, we would 
like to demonstrate it explicitly at the level of classical mechanics. 

Following Pars [63] , the motion of a point of unit mass in the field of gravity 
is described by the following equation 

f2 = {2Er^ + 2'yMr - a^)/r'^, (39) 

where a is the angular momentum integral (e.g. see Eq.(5.2.55) of Pars book). 
We would like now to replace r{t) by r{9) in such a way that dt = u{r{9))d9 
. Let therefore r{6) = exp(.T(^^)), -oo < x < oo. Unless otherwise specified, 
we shall write Tq = 1. In such (astronomical) system of units) we obtain, 

r = x' — exp(a;(0)). This result can be further simplified by choosing — 
exp(— a;(0)). With this choice Eq.(39) acquires the following form: 

{x' f = 2E + 27 Af exp(-a;) - exp(-2a;). (40) 

Consider points of equilibria for the potential U{r) ~ ~2jMr~^ + a'^r~'^ . From 

here we obtain: r* = According to Goldstein et al [80] such defined r* 

^^Not to be confused with the speed of hght ! 
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coincides with the major cUiptic semiaxis. It can be also shown, e.g. Pars, 
Eq.(5.4.14), that for the Kepler problem the following relation holds: E = 

. Accordingly, r* = , and, furthermore, using condition ^ = 

we obtam, — — = — or, a = — — . bmce m the chosen system ot 

units r(8) = ex.p{x{6)), we obtain, — ^ = exp(a;*(^)). It is conyenient to choose 

x*{6) = 0. This requirement makes the point x*{6) = as the origin and implies 
that with respect to such chosen origin = jM^^. Using this fact Eq.(40) can 
then be conveniently rewritten as 

^{x'f - 7M(exp(-a;) - ^ exp(-2a;)) = E (41a) 



or, equivalently, as 



y + A(exp(-2a;) - 2 exp(-a;)) = E, (41b) 



7M 

where A = Since this result is exact classical analog of the quantum 

Morse potential problem, transition to quantum mechanics can be done straight- 
forwardly at this stage. By doing so we haye to replace the Planck's constant 
h by clnA. After that, we can write the answer for spectrum at once [97]: 

- c\nA. 1,,. 

-^" = V[^-VW^"+2)] ■ 

This result contains an unknown parameter c which we would like to determine 
now. To do so it is sufficient to expand the potential in Eq.(41b) and to keep 
terms up to quadratic. Such a procedure produces the anticipated harmonic 
oscillator result ^ 

?^+Ax^ = E (43) 

with the spectrum giyen by En = (n+ ^)c\/2]41nv4. In the astronomical system 
of units the spectrum reads: En = {n+ ^)c2nhiA . This result is in agreement 
with Eq.(38). To proceed, we notice that in Eq.(38) the actual sign of the 
Planck-type constant is undetermined. Specifically, in our case (up to a con- 
stant) the energy En is determined by In^^^ = — In A so that it makes sense 
to write —En ~ nln A. To relate the classical energy defined by the Kepler-type 



^^In doing so some caution should be exercised since upon quantization equation r* = 

7M 

becomes r* = — ^ . Selecting the astronomical scale r J = 1 as the unit of length implies then 
that we can write the angular momentum as x — and to define >c as a| = a^. 
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7M 

equation E = to the energy we just have obtained, we have to replace this 

2r* 

Kepler-type equation by —En = — ln|i?[ = — 21n\/27r + lnr„ This is done in 
view of the 3rd Kepler's law and the fact that the new coordinate x is related to 
the old coordinate r via r = e^. Using Eq.(37) (for n = 1) in the previous equa- 
tion and comparing it with already obtained spectrum of harmonic oscillator 
we obtain: 

-2 In \/27r + In roB = -c2ir In A, (44) 

where in arriving at this result we had subtracted the nonphysical ground state 
energy. Thus, we obtain: 

1 27r2 

c=— !-^ln^. (45) 
27rlnA roB ^ ' 

Substitution of this result back into Eq.(42) produces 

-fi. = MI - (^;^)|- 2,ni - 5i^(n + i)P 

~27r^-4(n+i)+0.2(n+i)2. (46) 

To determine the number of bound states, we follow the same procedure as 
developed in chemistry for the Morse potential. For this purpose'*^ we intro- 
duce the energy difference A£^„ = — En = 4 — 0.4(n + 1) first. Next, 
the maximum number of bound states is determined by requiring Ai5„ = 0. In 
our case, we obtain: rimax = 9. This number is in perfect accord with observ- 
able data for planets of our Solar System (with Pluto being excluded and the 
asteroid belt included). In spite of such a good accord, some caution must be 
exercised while analyzing the obtained result. Should we not insist on physical 
grounds that the discrete spectrum must contain only finite number of levels, 
the obtained spectrum for the harmonic oscillator would be sufficient (that is 
to say, that the validity of the T-B law would be confirmed). Formally, it solves 
the quantization problem completely in accord with numerical data [93]. The 
problem lies however in the fact that these data were fitted to the power law, 
Eq.(37), in accord with the original T-B empirical guess. Heisenberg's honey- 
comb rule, Eq.(5), does not require the specific n— dependence. In fact, we 
have to consider the observed (the Titius-Bode-type) n— dependence only as a 
hint. With theoretical guidance emerging from this work, it is hoped, that the 
attempts will be made to fit the observational data to the Morsc-likc spectra 
in a way it is done routinely in chemical physics for the Morse-type potentials. 
In this work we intentionally avoid use of any adjustable parameters since the 
developed procedure, when supplied with correctly interpreted numerical data, 
should be sufficient for obtaining results without any adjustable parameters. 

Having said this, we must notice that there is still room for improving results 
we just obtained. Indeed, the constant c was determined using the harmonic 



Recall that in chemistry the Morse potential is being routinely used for description of the 
vibrational spectra of diatomic molecules. 
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approximation for the Morsc-typc potential. This approximation might fail very 
quickly as the following arguments indicate. Although we can calculate r'^s us- 
ing the T-B like law, Eq.(37), the arguments following this equation cause us to 
use the equation —En = —\n\E\ = — 21n\/27r + lnr„ for this purpose. This 
means that we have to use Eq.(46) (with ground state energy subtracted) in this 
equation in order to obtain the result for r„.If we ignore the quadratic correction 
in (46) (which is equivalent of calculating the constant c using harmonic oscilla- 
tor approximation to the Morse potential) then, by construction, we recover the 
T-B result, Eq.(46). If, however, we do not resort to such an approximation, 
calculations will become much more elaborate and arc not physically illuminat- 
ing. This is so because the T-B law, Eq.(37), is a purely empirical best fit to 
the observed data. In view of our calculations, Eq.(37) should be replaced by a 
more elaborate fitting result which is in agreement with data for the Morse-type 
potential. Since corrections to the harmonic oscillator potential in the case of 
the Morse potential are typically small, they do not change things qualitatively. 
Hence, we do not account for these complications in our paper. Nevertheless, 
accounting for these (anharmonic) corrections readily explains why the empir- 
ical T-B law works well for small n's and becomes increasingly unreliable for 
larger n's [91]. 

In support of our conjectures we performed similar calculations for satellite 
systems of Jupiter, Saturn, Uranus and Neptune. To do such calculations the 

astronomical system of units is not immediately useful since in the case of heavy 
planets one cannot use the relation Att'^/^Mq = 1. This is so because we have 
to replace the mass of the Sun Mq by the mass of the respective heavy planet. 

To do so, we write 47r^ = 7M0, multiply both sides by Mj (where j stands 
for the j-th heavy planet) and divide both sides by Mq. Thus, we obtain: 

4:TT^qj = jMj, where qj = —j- . Since the number is of order 10^^ — 10~^, 

it is inconvenient in actual calculations. To by pass this difficulty, we need to 
readjust Eq.(40) by rescaling x coordinate as a; = (5a; and, by choosing = qj. 
After transition to quantum mechanics such a rescaling results in replacing 
Eq.(42) for the spectrum by the following result: 

Since the constant c is undetermined initially, we can replace it by c = cS so 
that we reobtain back equation almost identical to Eq.(46). That is 

In this equation jMj = ^ir'^qj and {rj)i is the semimajor axis of the satellite 
lying in the equatorial plane and closest to the j-th planet. Our calculations are 
summarized in the Table 4 below. Appendix C contains the input data used in 
our calculations of n^^g^^^. 
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Table 4 



Satellite system\nniax 


n * 

theory 


Khs 


Solar system 


9 


9 


Jupiter system 


11-12 


8 


Saturn system 


20 


20 


Uranus system 


40 


18 


Neptune system 


33 


6 



Since the discrepancies for Uranus and Neptune systems may be genuine or not 
we come up with the following general pattern described below. 

4.4 Further analogies with atomic mechanics 

From atomic mechanics wc know that the approximation of independent elec- 
trons used by Bohr fails rather quickly with increased number of electrons. For 
this reason to expect that the T-B law is going to hold for satellites of all 
heavy planets is rather naive as wo explained already. At the same time, for 
planets rotating around the Sun such an approximation is seemingly good. The 
S0(2,l) symmetry explains why the motion of all planets should be planar but 
it does not explain why the motion of all planets is taking place in the plane 
almost coinciding with the equatorial plane of the Sun. The same is true for 
the regular satellites of all heavy planets as discussed by Dermott [95]. If we 
adopt the quantum mechanical point of view, then wc should accept that such 
an arrangement of planets is the result of some kind of spin-orbital interac- 
tion whose exact quantum mechanical nature remains to be elucidated. Other 
rotational resonances ubiquitous in the Solar system could then be explained 
quantum mechanically as well. The equatorial plane in which planets (satel- 
lites) move can be considered as some kind of an orbital (in the atomic physics 
terminology). It is being filled in accordance with the equivalent of the Pauli 
principle: eacJi orbit can be occupied by no more than one planet^^. Once the 
orbital is filled, other orbitals associated with other planes will begin to be filled 
out^''. Some of orbitals can be empty. This is indeed observed [95]. It should be 
said though that it appears (according to available data, e.g. see [92], that not 
all of observed satellites are moving on stable orbits. It appears also as if the 
"inner shell", when completely filled, acts as some kind of an s-type spherical 
orbital since the orbits of other sateUites he strictly outside the sphere whose 
diameter is greater or equal to that corresponding to the last allowed energy 
level in the Brst shell. The location of secondary planes appears to be quite 

*^The meteorite belt can be looked upon as some kind of a ring. We shall briefly discuss 
the rings below. 

^^Incidentally, such a requirement automatically excludes Pluto from the status of a planet. 
Indeed, although the T-B-type law, Eq.s(36),(37), can be seemingly adjusted to accomodate 
Pluto. Not only this would contradict the data summarized in Table 4 but also would be in 
contradiction with observational astronomical data for Pluto. According to these data the 
orbit inclination for Pluto is 17° as compared to the rest of planets whose inclination is within 
boundary margins of ±2°( except for Mercury for which it is 7°). 
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arbitrary as well as the filling of their stable orbits. Furthermore, without ac- 
count of spin-orbital interactions, quantum mechanics says nothing about the 
direction of orbital rotation. Although for all planets it does coincide with the 
direction of rotation of Sim's axis, in the case of Phoebe- the irregular satellite 
of Saturn-rotation takes place in the opposite direction to that of the axis of 
Saturn. If the spin-orbital interaction does exists, most likely, Phoebe's orbit is 
not a stable one. 

It is tempting to extend the picture just sketched beyond the scope of our 
Solar System. If for a moment we would ignore relativistic effects (they will be 
discussed in the next section) , we can then find out that our Sun is moving along 
almost circular orbit around our galaxy center with the period T = 185 • 10® 
years [98]. Our galaxy is also flat as our Solar System and the major mass is 
concentrated in the galaxy center. Hence, again, if we believe that stable stellar 
motion is taking place along the geodesies in accordance with laws of Einstein's 
general relativity, then we have to accept that our galaxy is a quantum object. 
It would be very interesting to estimate the number of allowed energy levels for 
our galaxy and to check if the Pauli-like principle works for the galaxy as well. 

4.5 Latest developments supporting our point of view 

We begin with the following observation. The motion of a planet of mass toq in 
the field of two static centers of attraction with masses mi and TO2 was discussed 
by Legendre and Jaeobi in 19th century [63] in connection with their study of 
elliptic functions. Such an idealized problem is a precursor of the restricted 
3-body problem to be discussed in the next subsection in connection with dy- 
namics of planetary rings. As simple as it is, the full study of this problem is 
extremely complex. It involves classification of all points and lines of equilibria 
and motions in the domains restricted by these lines. In addition to the eight 
major types of bounded orbits there are many more coming from collision of 
equilibrium point/lines etc. Characterization of the unbound motion is also 
interesting but is less complex. In quantum mechanics the motion of an elec- 
tron in the presence of two fixed positive ions is also a benchmark problem (in 
addition to study of He discussed in Section 1). Normally, charges of ions are 
assumed to be the same (e.g. for H^ ) which makes such a problem somewhat 
different (since they repel each other) from the problem studied by Legendre 
and Jacobi. All classification of molecular spectra can be traced back to this 
problem [97] . As in the case of H atom, the correspondence principle is not well 
established in this case since (to our knowledge) nobody studied the agreement 
between the quantum -mechanical calculations in the semiclassical limit and the 
results of Legendre-Jacobi theory modified due to the chemical requirements. 
Interestingly enough such a comparison was made to a larger extent between 
the classical restricted 3-body problem and its quantum analog. The quantum 
analog of the restricted 3-body problem exists in the form of the H atom placed 
in a strong crossed constant electric and magnetic fields [99]. Since semiclassical 
and classical analysis of such a system is sufficiently well understood, this fact 
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allows such a system to be studied both theoretically and experimentally. These 
studies are well summarized in two recent reviews [2,3] to which we refer our 
readers for details. For immediate purposes of this work the following quota- 
tion from Porter and Cvitanovic is helpful: " almost perfect parallel between 
the governing equations of atomic physics and celestial mechanics implies that 
the transport mechanism for these two situations is virtually identical: on the 
celestial scale, transport takes a spacecraft from one Lagrange point to another 
until it reaches its desired destination''''. On the atomic scale, the same type 
of trajectory transports an electron initially trapped near the atom across the 
escape threshold (in chemical parlance, across a "transition state"), never to 
return. The orbits used to design space missions thus also determine the ioniza- 
tion rates of atoms and chemical reaction rates of molecules" . This statement 
is nicely illustrated in the paper by Jaffe ct al [100] in which it is demonstrated 
that the transition state theory developed initially in chemistry (to describe 
the rates of chemical reactions) is working actually better in celestial mechanics 
where the discrepancy between the chemical theory and numerical simulations 
(done for celestial mechanics transport problems) is less than 1%. It should be 
noted though that the calculations were done at the classical level only (that 
is for a very large quantum numbers). The current status of transition state 
theory at the quantum and classical levels in chemistry can be found in the 
recent book by Micha and Burghardt [101].. 

4.6 The restricted 3-body problem and planetary rings 

Although the literature on restricted 3-body problem is huge, we would like 
to discuss this problem from the point of view of its connection with general 
relativity and quantization of planetary orbits along the lines advicated in this 
paper. 

We begin with several remarks. First, the existence of ring systems for all 
heavy planets is well documented [92]. Second, these ring systems are inter- 
spersed with satellites of these heavy planets. Third, both rings and satellites 
lie in the respective equatorial planes so that the satellites move on stable orbits. 
From these observations it follows that: 

a) While each of heavy planets is moving along the geodesies around Sun, 
the respective satellites are moving along the geodesies around 

respective planets: 

b) The motion of these satellites is almost circular. 

The restricted 3-body problem can be formulated now as follows. Given 
that the rings are made of some kind of small objects whose masses can be 
neglected''^ as compared to masses of both the satellite(s) and the particular 
heavy planet. Following previously discussed ideas by Laplace, we can ignore 
mutual gravitational interaction between these objects. Under such conditions 
we end up with the restricted three-body problem of motion of a given piece 

*'*E.g. see also paper by Convay at al [4]. 

*^This approximation is known as Hill's problem/approximation in the restricted 3-body 
problem [18]. 
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of a ring (of zero mass) in the presence of two bodies of masses mi and m2 
respectively. To simplify matters, one usually assumes that the motion of these 
two masses takes place on a circular orbit with respect to their center of mass. 
Complications associated with the eccentricity of such a motion are discussed in 
the book by Szebehely [102] and can be taken into account if needed. They will 
be ignored nevertheless in our discussion since we shall assume that satellites 
of heavy planets move on geodesies so that the center of mass coincides with 
the position of a heavy planet anyway thus making our computational scheme 
compatible with Einsteinian relativity. By assuming that ring pieces are mass- 
less we also arc making their motion compatible with requirements of general 
relativity since whatever orbits they may have-these are geodesies. 

Thus far only the motion of satellites in the equatorial planes (of respective 
planets) was considered as stable (and, hence, quantizable) . The motion of 
ring pieces was not accounted thus far by these stable orbits. The task now 
lies in showing that satellites lying inside the respective rings of heavy planets 
are essential for stability of these rings motion and, hence, they are making it 
quantizable. 

For the sake of space, we would like only to provide a sketch of arguments 
leading to such a conclusion. Our task is greatly simplified by the fact that a 
very similar situation exists for 3-body system such as Moon, Earth and Sun. 
Dynamics of such a system was studied very thoroughly by Hill whose work 
played profoimd role in Poincare' studies of celestial mechanics [88]. Recently, 
Avron and Simon [103] adopted Hill's ideas in order to develop a formal quan- 
tum mechanical treatment of the Saturn rings. In this work we follow the 
original Hill's ideas concerning dynamics of the Earth-Moon-Sun system. We 
claim that, when these ideas are looked upon from the point of view of modern 
mathematics of exactly integrable systems, they enable us to describe not only 
the Earth-Moon-Sun system but also the dynamics of rings of heavy planets. 
These modern mathematical methods enable us to find a place for the Hill's 
theory within general quantization scheme discussed in previous sections. 

4.6.1 Basics of the Hill's equation 

To avoid repetitions, we refer our readers to the books of Pars [63] , Chebotarev 
[98] and Brouwer and Clcmcncc [104] for detailed and clear account of the 
restricted 3-body problem and Hill's contributions to Lunar theory. Here we 
only summarize the ideas behind Hill's ground breaking work. 

In a nutshell his method of studying Lunar problem can be considered as 
extremely sophisticated improvement of previously discussed Laplace and La- 
grange method. Unlike Laplace, Hill realized that both Sun and Earth are 
surrounded by the rings of influcncc^^. The same goes for all heavy planets. 
Each of these planets and each satellite of such a planet will have its own do- 
main of influence whose actual width is controlled by the Jacobi integral of 
motion. For the sake of argument, consider the Saturn as an example. It has 

*® Related to the so called Roche limit [92]. 
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Pan as its the innermost satellite. Both the Saturn and Pan have their re- 
spective domains of influence. Naturally, we have to look for the domain of 
influence for the Saturn. Within such a domain let us consider a hypothetical 
closed Kcplcr-likc trajectory. Stability of such a Lagrangian trajectory is de- 
scribed by the Hill equation**^. Since such an equation describes a wavy- type 
oscillations around the presumably stable trajectory, the parameters describing 
such a trajectory are used as an input (perhaps, with subsequent adjustment) 
in the Hill equation given by 

(fx 

— + (go + 2gi cos 2t + 2g2 cosAt + - ■ ■)x = 0. (49) 

If we would ignore all terms except go first, we would naively obtain: XQ{t) = 
A(jcos{ts/qo+ e). This result describes oscillations around the equilibrium posi- 
tion along the trajectory with the constant go carrying information about this 
trajectory and the amplitude A is expected to be larger or equal to the average 
distance between the pieces of the ring. This naive picture gets very complicated 
at once should we use the obtained result as an input into Eq.(49). In this case 
the following equation is obtained 

(Px 

-^+<}oX + Aoqi{cos[t{y/q^ +2)+e] + cos[t(^ - 2) - e]} = (50) 

whose solution will enable us to determine grand Ai using the appropriate 
boundary conditions. Unfortunately, since such a procedure should be repeated 
infinitely many times, it is obviously impractical. Hill was able to design much 
better method. Before discussing Hill's equation from the perspective of modern 
mathematics, it is useful to recall the very basic classical facts about this 
equation summarized in the book by Ince [105]. For this purpose, we shall 
assume that the solution of Eq.(49) can be presented in the form 

oo 

x{t) = e"* E bre'''*- (51) 

r= — OO 

Substitution of this result into Eq.(49) leads to the following infinite system of 
linear equations 

oo 

(a + 2ri)% + J2 <lkK-k = 0, r e Z. (52) 

fe= — oo 

As in finite case, obtaining of nontrivial solution requires the infinite deter- 
minant A (a) to be equal to zero. This problem can be looked upon from 
two directions: either all constants qk are assigned and one is interested in the 
bounded solution of Eq.(51) for t — > oo, or one is interested in relationship 

*''ln fact, there will be the system of Hill's equations in general [98]. This is so since 
the disturbance of trajectory is normally decomposed into that which is perpendicular and 
that which is parallel to the Kepler's trajectory at a given point. We shall avoid these 
complications in our work. 
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between constants made in such a way that a = 0. In the last case it is impor- 
tant to know wether there is one or more than one of such solutions available. 
Although answers can be found in the book by Magnus and Winkler [106], we 
follow McKcan and Mocrbckc [107], Trubowitz [108] and Moscr [109]. 

For this purpose, we need to bring our notations in accord with those used 
in these references. Thus, the Hill operator is defined now as Q{q) = —■^ + q{t) 
with periodic potential q{t) = q{t + 1). Eq.(49) can now be rewritten as 

Q{q)x = Xx. (53) 

This presentation makes sense since Qq in Eq.(49) plays a role of A in Eq.(53). 
Since this is the second order differential equation, it has formally 2 solutions. 
These solutions depend upon boundary conditions. For instance, for periodic 
solutions such that x{t) = x{t + 2) the "spectrum" of Eq.(53) is discrete and is 
given by 

— oo < Ao < Ai < A2 < A3 < A4 < • • • t +00. 
We wrote the word spectrum in quotation marks because of the following. 
Eq.(53) does have a normalizable solution only if A belongs to the (pre assigned) 
intervals (Aq, Ai), (A2, A3), (A2i, A2i+i), ... In such a case the eigenfunctions 
Xi are normalizable in the usual sense of quantum mechanics and form an or- 
thogonal set. The periodic solutions make sense for vertical displacement from 
the reference trajectory. For the horizontal displacement the boundary condi- 
tion should be chosen as x(0) = x(l) = 0. For such chosen boundary condition 
the discrete spectrum also exists but it lies exactly in the gaps between the 
intervals just described, i.e. Ai < /xi < A2 < A3 < /X2 < A4 < • • •. For such 
a spectrum there is also set of normalized mutually orthogonal eigenfunctions. 
Thus in both cases quantum mechanical description is assured. One can do 
much more however. In particular, Trubowitz [108] designed an explicit pro- 
cedure of recovering the potential q(t) from the /i— spectrum supplemented by 
information about normalization constants. 

The Hill's equation can be interpreted in terms of the auxiliary dynamical 
(Neumann) problem. Such an interpretation is very helpful for us since it allows 
us to include the quantum mechanics of Hill's equation into general formalism 
developed in Sections 2 and 3. 

4.6.2 Connection with the dynamical Neumann problem and the 
Korteweg -de Vries equation 

Before describing such connenctions, we would like to add few details to the 
results of previous subsection. First, the number of the pre assigned intervals is 
always finite. This means that, beginning with some pre assigned i, we would be 
left with A2i = A2i+iVf > iThese double eigenvalues do not have independent 
physical significance since they can be determined by the set of single eigenvalues 
(for which A2i ^ A2i+i) as demonstrated by Hochstadt [110]. Because of this, 
the potentials q{t) in the Hill's equation are called the finite gap potentials**. 

*®Since there is only finite number of gaps [Ai, A2], [A3, A4], ...where the spectrum is forbid- 
den. 
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Hence, physically, it is sufficient to discuss only such potentials which possess 
finite single spectrum. The auxiliary /x— spectrum is then determined by the gaps 
of the single spectrum as explained above. With this information in our hands, 
we arc ready to discuss the exactly solvable Neumann dynamical problem. It is 
the problem about dynamics of a particle moving on the n— dimensional sphere 
< ^, ^ >= +■ ■ ■ + = 1 under the influence of a quadratic potential 
(j){£,) =< C,A^ > . Equations of motion describing the motion on n— sphere 
are given by 

C = -A^ + u{0^ with u{0 = m- <ii>- (54) 

Without loss of generality, wc assume that the matrix A is already in the diag- 
onal form: A := diag{ai, q„). With such an assumption we can equivalently 
rewrite (54) in the following suggestive form 

(^-^+u{m)j^k = ak^k ; k = l,...,n. (55) 

Thus, in the case if we can prove that u{S,{t)) in (55) is the same as q{t) in (53), 
the connection between the Hill and Neumann's problems will be established. 
The proof is presented in Appendix D. It is different from that given in the 
lectures by Moscr [109] since it is more direct and much shorter. 

This proof brought us the unexpected connection with hydrodynamics through 
the static version of Korteweg-de Vrics equation. Attempts to describe the Sat- 
urnian rings using equations of hydrodynamic are described in the recent mono- 
graph by Esposito [111]. This time, however, wc can accomplish more using 
just obtained information. This is the subject of the next subsection. 

4.6.3 Connections with SO(2,l) group and the K-Z equations 

Following Kirillov [112], we introduce the commutator for the fields (operators) 
^ and 1] as follows: [^,77] = ^drj — r]d^. Using the KdV, Eq.(D.lO), let us con- 
sider 3 of its independent solutions: ^0,^-1 and ^1. All these solutions can be 
obtained from general result: = t''''^^ + 0{t'^), valid near zero. Consider now 
a commutator [^0)^i]- Straightforwardly, we obtain, [^oi^i] = Analogously, 
we obtain, [^0;C-i] = and, finally, [^i,C-i] = —2^0- According to Kir- 

illov, such a Lie algebra is isomorphic to that for the group SL{2, R). Vilcnkin 
[113] demonstrated that the group SL{2,R) is isomorphic to SU{1, 1). Indeed, 

z — z 

by means of transformation: w — , it is possible to transform the upper 

z + 1 

half plane (on which SL{2,R) acts) into the interior of unit circle on which 
SU{1, 1) acts. Since, according to Appendix B, the group SU{1, 1) is the con- 
nected component of SO{2, 1), the anticipated connection with SO{2, 1) group 
is established. 

In Appendix D we noticed connections between the Picard-Fuchs, Hill and 
Neumann-type equations. In a recent paper by Veselov et al [114] such a con- 
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ncction was developed much further resuhing in the K-Z type equations for 
Neumann- type dynamical systems. We refer our readers to the original litera- 
ture, especially to the well written lecture notes by Moser [109]. These notes 
as well and his notes in collaboration with Zehndcr [84] provide an excellent 
background for the whole circle of ideas relating Hill's equation to integrable 
models. 

5 Solar System at larger scales: de Sitter, anti 
-de Sitter and conformal symmetries compat- 
ible with orbital quantization 

The obtained results demonstrate a remarkable interplay between the Newto- 
nian and Einsteinian mechanics already at the scale of our Solar System. Since 
quantization of stable orbits described in this paper is possible only with use of 
the basic experimental facts assuring correctness of results of general relativity, 
it is only natural to reverse this statement and to say that tie correctness of 
general relativity is assured by the observed pattern of stable (quantum) orbits. 

Since quantum mechanics can be developed group-theoretically, the same 
should be true for relativity. Quoting Einstein, Infeld and Hoffmann [89]: "Ac- 
tually, the only equations of gravitation which follow without ambiguity from 
the fundamental assumptions of the general theory of relativity are the equa- 
tions for empty space^°, and it is important to know whether they alone are 
capable of determining the motion of bodies" . In this work we argue that this 
is certainly correct locally when the Lorentzian-type symmetry holds true. Now 
we would like to discuss how such locally Lorentzian space-time embeds into 
space-times of general relativity possessing larger symmetry groups^^. Since 
this topic is extremely large, we shall discuss only the most basic facts from the 
point of view of results obtained in this paper. 

To our knowledge, Dirac [116] was the first who recognized the role of space- 
time symmetry in quantum mechanics. In his paper he wrote: "The equations 
of atomic physics are usually formulated in terms of space-time of special rel- 
ativity. They then have to form a scheme which remains invariant under all 
transformations which carry the space-time over into itself. These transforma- 
tions consist of the Lorentz rotations about a point combined with arbitrary 
translations, and form a group.... Nearly all of more general spaces have only 

■^^E.g. sec Eq.(29) of Section 2. 
^'^See also Section 6. 

^^At the level of quantum field theory Utiyama [115] demonstrated that the requirement 
of the local gauge invariance implemented for the non Abelian Lorentz group produces the 
Einstein field equations for gravitational field. This result implies that any "improvements" 
of Einsteinian relativity should involve changes in the local Lorentzian structure of space-time 
which is very unlikely. Independent arguments supporting this point of view are presented in 
Section 6. 
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trivial groups'''^of operations which carry the spaces into themselves.... There 
is one exception, however, namely the de Sitter space (with no local gravita- 
tional fields). This space is associated with a very interesting group, and so the 
study of the equations of atomic physics in this space is of special interest, from 
mathematical point of view." Subsequent studies indicated that the symmetry 
of space-time could be important even at the atomic scale [117,118]. Another 
reason to look at larger symmetry groups is associated with the cosmological 
constant problem [119] and, the associated with it problem of existence of 
cold dark energy (CDE) [120] cold dark matter (CDM) [121] and the modified 
Newtonian dynamics (MOND) [122]. Clearly, wc are unable to discuss these 
issues within the scope of this paper since they are more relevant to processes 
at galactic scales. Nevertheless, we would like to notice that, for instance, the 
MOND presupposes use of Newtonian and the modified Newtonian mechanics 
at the galactic scales which, as discussed in Section 4, strictly speaking, is not 
permissible even at the scales of our Solar System. The rationale for the dark 
energy and dark matter is explained in our recent paper [123] based on mathe- 
matical arguments consistent with that used by Grigory Perelman in his proof 
of the Poincare' conjecture. 

Hence, wc proceed with description of the de Sitter and anti-dc Sitter spaces 
based on results of our recent work. We begin with the following Hilbert- 
Einstein functional 

s'i9) = iM^^xRVg + ^iM^^^Vg (56) 

defined for some (pseudo) Riemannian manifold M of total space-time dimen- 
sion d. The (cosmological) constant A is determined as follows. 

Using Rij, the Ricci curvature tensor, the Einstein space is defined as solu- 
tion of the following vacuum Einstein equation 

Rij = Xgij (57) 

with A being a constant. Prom this definition it follows that 

R = dX. (58) 

At the same time, variation of the action S'^{g) produces 

Gij + ^Agij = 0, (59) 

where the Einstein tensor Gij is defined as Gij = Rij — \gijR with R being 
the scalar curvature determined by the metric tensor Qij . Combined use of 

®^This statement of Dirac is not correct. However, it is correct at the time of writing of his 
paper. 

^^Eq.(59) illustrates the meaning of the term "dark matter". The constant A enters into 
the stress-energy tensor (in the present case given by — ^Kgij) typically associated with the 
matter, Einstein [82]. 
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Eq.s(58) and (59) produces: A = X{d— 2). Substitution of this result back into 
Eq.(59) produces: 

,1 _ 1, 

Since by design Gj, ft* = 0, we obtain our major result: 



G} = h-:T)'5]^- (60) 



implying that scalar curvature R is constant. 

For isotropic homogenous spaces the Riemann curvature tensor can be pre- 
sented in the following known form [81]: 

Riiki = k{x){gikgji - guQjk)- (62) 

Accordingly, the Ricci tensor is obtained as: Rij = k{x)gij{d — 1). The Schur's 
theorem [124] guarantees that for d > 3 wc must have k{x) = k = const for the 
entire space. Therefore, we obtain: A = {d—l)k and, furthermore, R = d{d—l)k. 
The spatial coordinates can always be rescaled so that i? = fc or, alternatively, 
the constant k can be normalized to unity. For k > 0, k = and A: < 
we obtain respectively de Sitter, flat and anti-de Sitter spaces. Thus, we just 
have demonstrated that homogeneity and isotropy of space-time is synonymous 
with spaces being dc Sitter, flat and anti-dc Sitter very much like in ordinary 
Riemannian geometry there are spaces of positive, negative and zero curvature. 
This fact can be used to give the alternative description of just obtained results. 

We begin with simple observation that the surface of constant positive 
curvature is conformally equivalent to a sphere embedded in the Euclidean space 
[123] . In particular, let us consider a 3-sphere embedded into 4d Euclidean space. 
It is described by the equation 

= {xG Ei, xl + xl + xl + xl = R^}. (63) 

5^ is homogenous isotropic space with positive scalar curvature whose value 
is 6/R^. The group of motions associated with this homogenous space is the 
rotation group 150(4). The space of constant negative curvature is obtained 
analogously. For this purpose it is sufficient, following Dirac [116], to make 
xi purely imaginary and to replace R^ by — i?^ in Eq.(63). Such replacements 
produce: 

H^ = {xe Mi, x\-xl-x\-x\ = R^}. (64) 

In writing this result we have replaced the Euclidean space by the Minkowski 
space M4 so that the rotation group 5*0(4) is now replaced by the Lorentz group 
50(3, 1). The de Sitter space can now be obtained according to Dirac (1935) 
as follows. In Eq.(63) we replace E^ by E^, and make x\ purely imaginary thus 
converting into M5. The obtained space is the de Sitter space whose group 
of symmetry is 50(4, 1) 

dSi = {x e M5, x\- x\- x\- x\- x\ = R"^}. (65) 
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It has a constant positive scalar curvature whose value is 12/ B?. Very nice de- 
scription of such a space is contained in the book by Hawking and Ellis [125]. 
The connection between parameter R and the cosmological constant A is given 

by -R = sj^- The anti-de Sitter space is determined analogously as also dis- 
cussed by Hawking and Ellis and by Dirac. Specifically, it is given by 

adSi = {a: G £^3,2, x\ ~ xl - xl - xl + xl = i?^}, (66) 

where the five dimensional space £3,2 is constructed by adding the time-like 
direction to M4. Hence, the symmetry group of adSi is SO(3,2). All these 
groups can be described simultaneously if, following Dirac [116], we introduce 
the quadratic form 

j:x^x^ = B? (67) 

in which some of the arguments are allowed to be purely imaginary. Trans- 
formations preserving such a quadratic form are appropriate respectively for 
groups S0(5), S0(4,l) and SO(3,2). We still can embed all these groups into 
a larger (conformal) group SO(4,2) by increasing summation from 5 to 6 in 
Eq.(67). In such a case all groups discussed in this work, starting from S0(2,l), 
can be embedded into this conformal group as subgroups as discussed in great 
detail by Wybournc [70]''''*. Comprehensive group-theoretic description of the 
Einstein spaces, e.g. see Eq.(57), including those which are invariant with re- 
spect to the conformal group, can be found in the monograph by Petrov [127]. 
The significance and use of c;onformal symmetry in both gravity and confor- 
mal field theories has been recently extended in [123]. All existing cosmological 
models in the limit 00 approach one of the Einstein's spaces whose group 
of symmetry belongs to the types just described. The de Sitter and anti-de 
Sitter spaces are the simplest examples of such spaces [128].^^. 

The task still remains to find out if representations of these larger groups, e.g. 
see Vilcnkin [113] and Wybourne [70] for mathematical details, can produce 
the exact solutions of radial Schrodinger equations not listed in the Natanzon- 
style classification, e.g. see Levai [72], for S0(2,l). If such solutions do exist, one 
might be able to find those of them which are of relevance to celestial quantum 
mechanics and, hence, to cosmology. 



6 Conclusions 

It is a remarkable historical fact (discussed in Section 4.2.) that Laplace was 
the first who studied resonance dynamics of known satellites of Jupiter (ef- 
fectively) using geodesies while Lagrange analyzing motion of these satellites 

''"^Incidentally, the work by Graner and DubruUe [126], when translated into group-theoretic 
language, becomes just a corollary of conformal invariaiice implied by the conformal group . 

■''■''The most recent mathematically rigorous description of both do Sitter and anti-de Sitter 
spaces can be found in the paper by Andersson and collaborators [129]. 
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along stable geodesies had arrived (effectively) at the Bohr-Sommerfeld quanti- 
zation condition. The validity of the geodesic-type approximation is based on 
the equivalence principle of general relativity correct in the limit of vanishingly 
small masses as compared to the mass of the central body. Using this principle 
Einstein [82] was able to calculate the perihelium shift for Mercury. Within ap- 
proximations he made all interactions of Mercury with the rest of planets were 
ignored. Because of this, the same type of calculations in the spirit of Laplace 
can be made for all planets as discussed by Misner et al [83]. The fundamental 
problem then lies in proving that such stable motions will survive in the case 
if masses of planets are small but nonzero. In the case of satellites of Jupiter 
such a task was completed by Poincare' and de Sitter as mentioned in Section 
4. The latest advancements are also discussed in the same section. 

In this (concluding) section we argue that it is possible to arrive at the 
same field equations of general relativity by entirely by passing the equivalence 
principle. This can be achieved by studying the limiting case of dynamics of 
2-f-l gravity as discussed in our papers [130-132] and in more recent paper by 
't Hooft [133]. In such a limit one studies surface dynamics of the fictitious 
2 dimensional gravitating bodies^^. Simple topological arguments applied to 
this case indicate that the Einstein field equations survive such a reduction 
while Newtionian (actually, the Poissonian-type) equations do not survive this 
reduction. As result, the dynamics of 2+1 gravity is strictly Einsteinian. In 
mathematics this type of dynamics ( the dynamics of measured foliations) was 
discovered totally independently of gravity-related considerations by Thurston 
[134] in his study of 3-manifolds^^. Physically, such type of dynamics is realized 
in dynamics of some 2 dimensional liquid crystals. All this is explained in our 
works [130,131] which, in turn, were inspired by the earlier work by Deser, 
Jackiw and t'Hooft [135]. The theory of foliations (for surfaces) is thoroughly 
discussed in monograph by Nikolaev [136]. The book by Hehl and Obukhov [137] 
uses foliations for description of classical electrodynamics in 3-1-1 space. Some 
basics of foliation theory from the point of view of Lie groups and Lie algebras are 
discussed in the very readable book by Moerijk and Mrcim [138]. By reversing 
reduction arguments it is possible to arrive at dynamics of full 3-1-1 gravity in the 
presence of matter beginning from the dynamics of measured foliations for 2+1 
gravity. The topological properties of dynamics of 2+1 gravity are described 
in terms of polynomials of knots and links as explained in detail in the paper 
by Kholodenko [132]. Accordingly, dynamics of full 3+1 gravity should be 
associated with time dynamics of 3-manifolds foliating 3+1 space [123]. These 
remarks provide needed physical justification to works by Witten [138] and 
Kholodenko [139] connecting statics and dynamics of 3 (or 2+1) gravity with 
conformal and string theories. 

Acknowledgements The Tables 1 and 2 are reproduced with permis- 

®^ These can be actually visualized by the ciosscctions (lying on the surface) of infinitely 
long and thin massive 3 dimensional rods. According to the imposed rules we are allowed to 
watch only the motion of crossections. 

^^He was awarded the Fields medal for these studies. 
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Appendix A. Details of Heisenberg's derivation of the commutator 
identity [x,p] = ih 

In this appendix we would like to provide some details of Heisenberg's rea- 
soning leading to the discovery of [x, p] = ih. This would be unnecessary should 
his original paper [27] contain all details. 

At the classical level consider a gas of noninteracting atoms, better just one 
atom containing N electrons which are assumed to scatter light independently. 
The interaction between the incoming light and individual electron is described 
with help of the combination d= /3E where d is the dipole moment of the 
electron in the atom, E is the strength of the external electric field which is 
assumed to be time-dependent, and /? is the polarization tensor (in the simplest 
case it is assumed to be a scalar). In the medium the strength of the electric 
field changes as compared to the vacuum. By denoting it as D it is known that 
D=E -|-47rP where P=Nd. Since, at the same time, by definition, d=er we 
have to have an equation for r. It is given by 

r + uj^r + -fi- =—E{t) (A.l) 
m 

where e is electron's charge and m is its mass. In writing this equation it is 

assumed that our electron is bound harmonically (with the basic frequency loq) 
and that the friction is of known (electromagnetic) nature and is assumed to 
be small. Using Fourier decomposition of r(t) we obtain, 

rH = -^ — Ix^- M 

This equation allows us to obtain P and, hence, D as follows : 

1 

D = E -I- 47rP = (1 47riV 5 5 )E = £(w)E. (A.3) 

This equation defines a complex frequency-dependent dielectric constant s{u). 
From electrodynamics it can be equivalently rewritten as e{Lu) = (n{uj)^iK{uj))'^ 
where n(a;) is the refractive index while >c(oj) is the coefficient of absorption. 
Using these facts we can write approximately 



n{u;) = l + 2nN --^—^y—. (A.4) 



1 

muiQ — u}'^ + iivy 

By ignoring friction in the high frequency limit we obtain, 

„2 



n{oj) = l-2wN (A.5.) 
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To account for quantum mechanical effects, Tliomas, Reich and Kuhn in 1925 
(just before the quantum mechanics was born !) have suggested to replace 
Eq.(A.4) by 

nM = l + 2.7V^5:-^ (A.6) 



where, following these authors, we ignored friction and introduced the oscillator 
strength /i.To reconcile Eq.(A.6) with (A. 5) we have to require X^/i = 1. This 

i 

requirement is known as the sum rule. These facts were known to Kramers and 
Heisenberg^^ where our readers can find additional details. To make our point 
and to save space, we would like to reobtain the result, Eq.(A.6), quantum 
mechanically using modern formalism. We refer our reader to the book by 
Davydov [140] for additional details. Basically, we need to calculate quantum 
mechanically the dipole moment d, that is 



eriprnd^r. (A.7.) 



In this expression the wave function Tpm is calculated with help of the stationary 
perturbation theory with accuracy up to the first order in perturbation (which 
is er-E). A short calculation produces the following result for the oscillator 
strength, 

This result can be equivalently rewritten as 

fkm = |(fc I X I m)* {k \ X \ m) + {k \ X \ m)* {k \ x \ m)}. (A. 9) 

Since, however, 

imujkm {k\x\m) = {k\px\m) (A. 10) 

we can rewrite Eq.(A.9) as 

fkm = \x\k){k\p:^\m) - {m\px\k) {k\x \ m)} (A.ll) 

in 

since ujkm = —^mk- Finally, we have to require J2fkm = 1- This is possible 

k 

only if 

— (m I xpx - PxX I m) = 1, (A.12) 
in 

QED. 

Appendix B. Some quantum mechanical problems associated with 

the Lie algebra of SO(2,l) group 



®*E.g. see the reference in Heisenberg's paper. 
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Following Wybourne [70] let us consider the second order differential equa- 
tion of the type 

fL + vix)Yix)=0 (B.l) 

where V{x) = a/x'^ + bx'^ + c. Consider as well the Lie algebra of the noncompact 
group S0(2,l) or, better, its connected component SU(1,1). It is given by the 
following commutation relations 

[Xi, X2] = -iXs; [X2, X3] = iX,; [X3, X,] = iX^ (B.2) 

We shall seek the realization of this Lie algebra in terms of the following gener- 
ators 

d d^ 
X^ := — + ai{x)] :=i[fc(a;)— +02(x)]; X3 := — + 03(2;). (B.3) 

The unknown functions ai{x) , a2{x) ^ aj,{x) and k{x) arc determined upon sub- 
stitution of Eq.s(B.3) into Eq.s(B.2). After some calculations, the following 
result is obtained 

d^ a x^ -i, d 1, a x"^ ,^ ^. 

^^^=d^ + ^+16^^^^=T[^^ + 2]^^^^=^ + ^-16- (^-^^ 

In view of this, Eq.(B.l) can be rewritten as follows 

[(^ + 86)Xi + (i - 86)X3 + c]Y{x) = (B.5) 

This expression can be further simplified by the unitary transformationC/XiC/~^ = 
Xi cosh 61 + X3 sinh6'; UXsU-'^ = Xi smhO + X3 coshe* with U = expi-iOX^). 
By choosing tanh6l = -(1/2 + 86)/(l/2 - 86) Eq.(B.5) is reduced to 

X:,Y{x) = -^Y{x) (B.6) 

where the cigenfunction Y{x) = UY{x) is an eigenfunction of both X3 and the 
Casimir operator = X| — X| — X^ so that by analogy with the Lie algebra 
of the angular momentum we obtain, 

X2Y,„(x) = J(J+ l)fj„(x) and (B.7a) 
XsYjnix) = -^Yjn{x) = i-J + n)Yj4x); n = 0, 1, 2, .... (B7b) 

It can be shown that J(J -|- 1) = —a/4 — 3/16. Prom here we obtain : J = 
— i(l ± ^ J — a); -J — a > 0. In the case of discrete spectrum one should choose 
the plus sign in the expression for J. Using this result in Eq.(B.7) we obtain 
the following result of major importance 



4n + 2 + Vl-4a= —=. (B.8) 
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Consider now the planar Kcplcr problem. In this case, in view of Eq.(32), the 
radial Schrodinger equation can be written in the following symbolic form 



Id V u 
r dr r 



9 



R{r) = 



(B.9) 



By writing r = and R{r) = x 2 TZ{x) This equation is reduced to the canon- 
ical form given by Eq(B.l). e.g. to 



, 4u + 1/4 ,0 , . 



(B.IO) 



so that the rest of arguments go through. Analogously, in the case of Morse-type 
potential we have the following Schrodinger-type equation initially: 



dz'^ 



■ pe 



k 



R{z) = 



(B.12) 



By choosing z = Inx^ and R{z) = x iTZ{x) Eq.(B12) is reduced to the canonical 
form 



Aa'^x^ 



By analogous manipulations one can reduce to the canonical form the radial 
equations for Hydrogen atom and for the 3-dimensional harmonic oscillator. 

Appendix C. Numerical data used for claculations of n-^heory 
(Table 4). 

1 au=149.598-10^fcm 

Masses (in kg): Sun l.QSS-lO^o, Jupiter 1.8986- lO^^, Saturn 5.6846-1026, 
Uranus 8.6832-1025, Neptune 10.243-10^^. 

qj : Jupiter 0.955-10-3, Saturn 2.86-10-^ Uranus 4.37-10-^ Neptune 5.15-10"^. 
(rj)^ (km) : Jupiter 127.69-10^, Saturn 133.58-10^, Uranus 49.77-10^ 
Neptune 48.23- 10^. 

( ) : Earth 4.0062, Jupiter 3.095, Saturn 1.844, Uranus 0.9513, 



In, 

V2ri 

Neptune 1.15. 



Appendix D. Connections between the Hill and Neumann's 
dynamical problems. 

We follow our paper [141] where some mathematical of the results of the pa- 
per by Lazutkin and Pankratova (1975) were used for solution of concrete phys- 
ical problems. In particular, following our paper, let us consider the Fuchsian- 
type equation given by 

y" + I'py = 0, (D.i) 
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where the potential cfi is determined by the equation cj) = [/] with / = 
and 1/1,2/2 being two independent solutions of Eq.(D.l) normalized by the re- 
quirement 2/^2/2 -2/2 yi = l.The symbol [/] denotes the Schwarzian derivative of 
/. Such a derivative is defined as follows 

f f" — ^ ( f'S^ 
[f\ = • (D-2) 

Consider Eq.(D.l) on the circle and consider some map of the circle given by 
F{t+1) = F{t) + l. Let t = F{£) so that yif) = Y{i)^F'{S) leaves Eq.(D.l) form 
-invariant, i.e. in the form Y" -\- — with potential $ being defined now as 
$(^) = (j){F{!^))[F'{£)Y + [-^"(0]- Consider next the infinitesimal transformation 
= 6 + with 5 being some small parameter and ip{C) being some 

function to be determined. Then, + 5ip{i,)) = (^(0 + S{fip){^) + 0{S^). 
Here (T'<^)(^) = <p{Ov'{^) + i</'"'(0 + 2<^'(0</='(0- Next, we assume that the 
parameter 6 plays the same role as time. Then, we obtain 

lim ^ = ^ = i<^'"(0 + m^'iO + '^'P'iOviO (D.3) 

Since thus far the perturbing function (p(^) was left undetermined, we can choose 
it now as (p{^) = ^(0- Then, we obtain the Korteweg -de Vriez (KdV) equation 

determining the potential ^{^). For reasons which will be explained in the text, 
it is sufficient to consider only the static case of KdV, i.e. 

r{o+m)^'{o = o. (D.5) 

We shall use this result as a reference for our main task of connecting the Hill 
and the Neumann's problems. Using Eq.(54) we write 

u{0=m-<ti>- (D.6) 

Consider an auxiliary functional (p{^) =< ^, A~^^ > . Suppose that <^(^) = w(^). 
Then, 

fjfl 

-=2<^,A^>-2<^,^> . (D.7) 

But < ^, ^ >= because of the normalization constraint < ^, ^ >= 1. Hence, 

^ = 2 < ^, A^ > . Consider as well ^ . By using Eq.s (54) it is straightforward 
dt dt 

to show that ~ ^ < ^, A~^^ > . Because by assumption ip{^) = u{^) we have 

to demand that < ^, A~^^ >=< ^, > as well. If this is the case, consider 
furthermore 

-^ = 2<^,A-'^>+2<^,A-'^> (D.8) 
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Using Eq.s(54) once again we obtain 



(fu 

= -2 + 2u(p + 2<^,A-^^> . (D.9) 

Finally, consider as well Using Eq.(D.9) as well as Eq.(54) and (D.7) we 
obtain, 

cPu ^du . du „ du ,^ 

= 2—ip + Au— = Qu— (D.IO) 
dt^ dV dt dt ^ ' 

By noticing that in Eq.(D.5) we can always make a rescaling (^(^) — > A<?f)(^) 
wc always can choose A = —1. Therefore Eq.s (D.5) and (D.IO) coincide. This 
establishes the correspondence between the Neumann and Hill-type problems. 
QED 
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